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SUMMARY 


The sonic boom flight test program conducted at Jackass Flats, Nevada, during the summer 
and fall of 1970 consisted of 121 sonic-boom-generating flights over the 466 m (1529-ft) BREN 
tower. This test program was designed to provide information on several aspects of sonic boom, 
including caustics produced by steady flight near the threshold Mach number, caustics produced by 
longitudinal accelerations, sonic boom characteristics near lateral cutoff, and the vertical extent of 
shock waves attached to near-sonic (M < 1 .0) airplanes. The measured test data (except for the 
near-sonic flight data) were analyzed in detail to determine sonic boom characteristics for these 
flight conditions and to determine the accuracy and the range of validity of linear sonic boom 
theory. 

The observations near caustics are of particular value since no methods are currently available 
for predicting realistic pressure signatures at caustics because of the nonlinear effects that 
predominate there. Overpressure increases measured at caustics produced during threshold Mach 
number flight, compared to the overpressure produced during steady, level flight at about Mach 1 .2, 
were relatively low, ranging from about 1.0 to 1.8, except in one case where it appeared that 
small-scale atmospheric turbulence produced “spikes” on a caustic signature. Caustics were also 
produced by small inadvertent changes in the airplane speed during several of the “steady,” level 
threshold Mach number flights. These caustics were slightly stronger, with a maximum amplification 
factor of about 3. Measured overpressures at caustics produced by airplane accelerations beginning 
from Mach 0.95 ranged from 2 to 5 times those which would be observed during steady, level flight 
at about Mach 1 .2. Pressure signatures were also observed near lateral cutoff which resembled those 
measured at caustics. These disturbances were of very low intensity, however, less than one-half the 
intensity beneath the flight path. The caustic phenomena are analyzed and documented in detail. 
The distinguishing features of pressure signatures near caustics are the “U” shape of the signature, 
about a 40% longer duration than normal, and the sharp peaks at the bow and tail shocks. 



Comparison of theoretical calculations with the observed data showed good agreement in all 
cases where it was possible to make such calculations. Shock wave intensities, shock wave 
inclination angles, and sonic boom arrival times were calculated and compared with the 
observations. The accuracy of the theoretical calculations was limited by the accuracy of the input 
airplane flight track and meteorological data, particularly for the low supersonic flights where the 
shock wave travels long distances almost parallel to the ground. Inaccuracies in the upper level wind 
conditions are also considered to be an important source of error. Caustic locations during transonic 
acceleration and lateral cutoff locations can be predicted to within ±1.0 km (3300 ft). Shock wave 
intensities agree reasonably well, as do signature shapes when the effects of small-scale turbulence 
are neglected in the observed data. Shock wave arrival times can be predicted to within ±1 .0 sec. 
Thus, the comparison of theory and experiment has tended to verify the theory and has also 
indicated its range of validity. The linear theory is invalid within about 100-200 m (330-660 ft) 
vertically above caustics and where the shock wave is within a few degrees of the cutoff condition. 

Analysis of the rumble data produced during flight near the threshold Mach number showed 
that “low rumbles” occurred when the airplane ground speed was at least 6 m/sec (19.7 ft/sec) 
lower than the maximum shock propagation speed. This, roughly, is the “safety factor” (or 
reduction in allowable airplane ground speed to avoid objectionable noise at the ground) that would 
be required for these flight and meteorological conditions. Comparison of a theoretical “safe 
altitude” for sonic boom cutoff with the observed data was good considering the assumptions made 
in deriving it. 


INTRODUCTION 


During the summer and fall of 1970 a series of sonic boom flight tests were conducted by the 
NASA at Jackass Flats, Nevada. During this flight test program, 121 sonic-boom-generating flights 
were made. These flights were designed to provide information on several aspects of sonic boom, 
including caustics produced by steady flight near the threshold Mach number, caustics produced by 
longitudinal accelerations, sonic boom characteristics near lateral cutoff, and the vertical extent of 
shock waves attached to near-sonic (M< 1.0) airplanes. By use of the 466 m (1529-ft) BREN 
tower, sonic boom signature measurements as a function of altitude were obtained for the first 
time. The primary goal of this test series was to obtain definitive data on caustics produced by 
accelerations and by atmospheric refraction (threshold Mach number and lateral cutoff). The 
analysis of these sonic boom data is the subject of this report. The test results and caustic 
phenomena are presented and described, and wherever possible comparisons are made with the 
linear sonic boom theory. 


The test arrangements for the threshold Mach number flights, the longitudinal accelerations, 
and the lateral cutoff flights are summarized in table 1 and schematically illustrated in figure 1. The 
desirabUity of using a tower to observe the sonic boom pressure field in the vertical plane is 
illustrated in figure 1, particularly for shock waves near cutoff. For a shock wave angle of incidence 
of 0.087 rad (5^, a ground array about 5.15 km (3.2 st mi) long would be required to obtain the 
same data sample as that measured on the 466 m (1529-ft) BREN tower; at 0.035 rad (2f), a linear 
array of 13.85 km (8.6 st mi) is required. In addition, the use of a tower gives a much better chance 
of observing the very localized caustic phenomena. The success of these tests in observing caustics 
can largely be attributed to the use of the BREN tower. 

During the last two decades a large amount of work has been devoted to the development of 
theoretical methods for predicting sonic boom characteristics. Fundamental contributions were 
made by Friedrichs, Hayes, Landau, and Whitham (refs. 1 to 6). The current state of the art in sonic 
boom prediction is represented by the comprehensive analysis and computer program assembled by 
Hayes, Haefeli, and Kulsrud (ref. 7). Although the basic theory is linear in nature (except for the 
calculation of pressure signature “aging” and shock locations), a wide range of effects can be 
accounted for, such as the effects of airplane configuration and lift, nonstandard meteorological 
conditions, and airplane maneuvers. This basic theory has been verified by wind tunnel tests and 
flight test data (ref. 8). 

During the last few years research efforts have been directed toward extending the basic theory 
to predict pressure levels in the vicinity of caustics, which can be produced by airplane maneuvers 
and/or by atmospheric refraction. A caustic is a point in space where ray focusing occurs and the 



TABLE 1.-SUMMARY OF 1970 BREN TOWER TESTS 


Type of flight 

Approximate 
altitude 
above ground, 
km (ft) 

Approximate 

Mach 

number 

Number 

of 

flights 

Steady, level flight 

9.14 

1 .07 to 1 .20 

79 

near the threshold 

(30 000) 



Mach number 




Longitudinal acceleration 

9.14 

0.95 to 1.3 

19 

(at constant altitude) 

(30 000) 



Lateral cut-off 

8.84 

1.3 

9 

(steady, level flight) 

(29 000) 



Near-sonic 

0.91 

0.95 to 1.00 

14 

(steady, level flight) 

( 3 000) 






121 


ray-tube area of linear sonic boom theory is zero. Linear theory predicts infinite sonic boom 
intensity at this point. The concept of a ray-tube area, however, is linear and does not account for 
the nonlinear effects that predominate close to caustics. 

The mathematical formulation of the pressure signature behavior near a caustic has been 
outlined by Guiraud (ref. 9) and Hayes (ref. 10). Th&y (ref. 11) has attempted to carry the 
mathematics to completion, but his results must be viewed as giving only a qualitative indication of 
the pressure amplification. Seebass (ref. 12) also attempted a mathematical solution for an acoustic 
disturbance (no shocks) entering a caustic region, while Coakley (ref. 13) has attempted a numerical 
solution. The calculation of sonic boom pressure signatures and shock intensities in a caustic region 
has not been accomplished at this time. However, flight test measurements have been most valuable 
in identifying the phenomena and will, in the future, serve to confirm theories that may be 
advanced. 

During the last decade there have been a number of flight test programs that have contributed 
substantially to the current understanding of sonic boom phenomena near caustics. One of the 
earliest was conducted at NASA, Wallops Station, where seven threshold Mach number flights were 
flown (ref. 14). The tests at Edwards AFB in 1961 consisted of several different airplane maneuvers, 
including four longitudinal accelerations (refs. 15-17). Another series at Edwards AFB in 1964 
included five threshold Mach number flights and five longitudinal accelerations (ref. 18). Other 
programs have provided information on the variability of sonic boom due to atmospheric effects as 
a function of lateral displacement from the flight track (refs. 19-23). Reference 24 contains a 




general discussion of the variability of sonic boom and some preliminary analysis of the Jackass 
Flats lateral cutoff data, and reference 25 contains a preliminary analysis of the Jackass Flats 
threshold Mach number flights. 

The French flight tests (ref. 26) demonstrated that the method of linear geometric acoustics is 
satisfactory for predicting caustic locations. These test results also demonstrated that a dense array 
of microphones was necessary for observing the caustic phenomena, since they occur over small 
ground areas. 

These early test series, while of excellent quality, did not give definitive measurements on 
certain aspects of sonic boom phenomena because of the lack of systematic measurements and the 
lack of the required density of microphones. One aspect that needed further investigation was the 
sonic boom characteristics near lateral cutoff, since several investigators have postulated large 
pressure magnifications at lateral cutoff (refs. 27-29). Another was the caustic and acoustic 
disturbances associated with threshold Mach number flight where a caustic (and cutofO is produced 
at some distance above the ground. Finally, the caustic produced by the transonic acceleration of 
supersonic airplanes is of particular interest since it will be produced by SSTs. These considerations 
prompted the 1970 test program at the BREN tower. 

The data measured by the NASA at the BREN tower represent a unique and valuable set of 
information on the some of the least well understood aspects of sonic boom propagation. Vertical 
measurements through the shock waves were made for the first time in any flight test program, and 
these resulted in the most definitive information available on the formation and nature of caustics. 
The data acquired during the test program are also valuable because they have application to several 
future airplane concepts and operations. Among them are: 

• Commercial transport operation near the threshold Mach number (M >1.0) to avoid 
sonic boom noise on the ground 

• Commercial transport operation at near-sonic (M < 1.0) speeds 

• SST operation during the transonic acceleration phase of flight 

• Pressure signature characteristics at the outer edge of the sonic boom carpet during 
normal SST operations 

In general, the objectives of the detailed analysis of these test data were: (1) to compare sonic 
boom theory with observations, and (2) on the basis of these comparisons, to interpret the 
measurements for each flight condition. The presentation of the test results and analysis in this 



report begins with the threshold Mach number data, followed by the longitudinal acceleration and 
the lateral cutoff data. (The transonic data obtained by flights at 900 m (3000 ft) altitude and Mach 
numbers less than 1 .0 were not analyzed in this study.) A brief discussion of the measurement 
procedure and the method used to calculate theoretical pressure signatures is given. The appendixes 
contain an evaluation of the effect of water vapor on shock propagation speed, tabulated airplane 
F-function data, a discussion of theoretical relationships near cutoff during threshold Mach number 
flight, and the derivation of a theoretical safe cutoff altitude during threshold Mach number flight. 


SYMBOLS 


This section defines symbols used in the analysis. 

Symbol 

A airplane acceleration 

Ajj ray-tube area projection onto horizontal plane 

A(S) ray-tube area normal to rays 
a sound speed 

a* sound speed using virtual temperature, Ty 

Cl hft coefficient 

Cq Snell’s law invariant 

D virtual temperature increment, (Ty - T); airplane lateral displacement from tower 

F airplane F-function for pressure signature calculations 

F^ geometry contribution to F-function 

F^j nondimensionalized geometry contribution 

Fg lift contribution to F-function 

Fgj nondimensionalized lift contribution 

Fj input F-function, sum of lift and geometry contributions 

Ff F-function conversion factor 

g acceleration of gravity 


h airplane altitude above mean sea level 

Kr ground reflection coefficient 

L distance from airplane nose 

Lp reference F-function length used for nondimensionalizing L 

M airplane Mach number = V/a 

Mp threshold Mach number 

i? unit vector normal to wave front 

np lift (normal) load factor 

nj thrust (axial) load factor 

P pressure 

q dynamic pressure 

r water vapor mixing ratio 

S separation distance between leading and trailing shocks produced by acceleration 
Srpf airplane reference wing area 

T air temperature 

tRpF reference time used for shock wave profile calculations 
virtual air temperature 
t^ shock wave arrival time 

u horizontal wind speed 

Uj^ component of horizontal wind in plane of the normal to the shock wave 
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Uf component of horizontal wind in plane tangent to the shock wave 

V airplane velocity relative to atmosphere, Ma^ 

Vq airplane ground speed, (Ma^ - Ujj^) 

Vp shock propagation velocity, (a - Uj^) 

V-, shock propagation speed determined from shock wave arrival times over ground microphone 
Pmic 

array 

V„„ shock propagation speed normal to the shock wave surface 

Vp* shock propagation velocity with increase because of water vapor, (a* - Uj^) 

W airplane weight 

(X,Y,Z) reference coordinate system: east, north, and above ground, respectively 

distance downtrack of acceleration caustic-ground intersection 

sonic boom cutoff altitude during threshold Mach number flight 

Zg Theoretical minimum altitude above the ground for which cutoff can occur to obtain 
low-intensity acoustic-like disturbances at ground during threshold Mach number flight 

a shock wave angle of incidence 

olq ground slope 

/3 Prandtl-Glauert parameter, (M^ - 1)^^^ 

A perturbation from undisturbed value 

€ ratio of molecular weight of water vapor to that of dry air 
1? direction from which wind blows 

0 inclination angle of wave normal n below horizontal 



/i Mach angle, sin"^ (1/M) 

V heading angle of wave normal n 

0 azimuth angle of wave normal from vertical plane 

0 airplane heading angle 

Subscripts 

G ground level 

max maximum 

REF reference value 

T tower 

o initial value at airplane altitude 
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GENERAL TEST DESCRIPTION 


The unique feature of the Jackass Flats sonic boom tests was the use of the 466-m (1529-ft) 
BREN tower shown in figure 2. Fifteen microphones were placed at 30.5 m (100-ft) intervals on the 
tower for observing the sonic boom pressure field. Fourteen microphones were placed on the 
ground in a line parallel to the nominal flight path. Two additional microphones were placed along a 
line normal to the flight path. Figure 3 schematically shows the test arrangement and table 2 gives 
the detailed microphone locations with respect to the tower base. The ground slopes slightly 
upward by about 0.017 rad (1“) from microphones G-1 to G-14. 


TABLE 2.-TOWER AND GROUND MICROPHONE LOCATIONS 



Tower microphones 

1 Ground microphones | 

Tower 

microphone 

number 

Altitude above tower base^ 

Ground 

microphone 

number 

Horizontal distance 
from tower base 

Elevation with 
respect to tower 
base 

m 

ft 

m 

ft 

m 

ft 

Tower Top 

466.10 

1529.2 

G-1 

-487.68 

-1600 

-10.27 

-33.7 

T-15 

458.36 

1503.8 

G-2 

-426.72 

-1400 

- 8.69 

-28.5 

T-14 

428.09 

1404.5 

G-3 

-365.76 

-1200 

- 7.32 

-24.0 

T-13 

397.40 

1303.8 

G^ 

-304.80 

-1000 

- 6.00 

-19.7 

T-12 

366.92 

1203.8 

G-5 

-243.84 

- 800 

- 4.97 

-16.3 

T-11 

336.47 

1103.9 

G-6 

-182.88 

- 600 

- 3.72 

12.2 

T-10 

305.96 

1003.8 

G-7 

-121.92 

- 400 

- 4.75 

-15.6 

T-9 

275-48 

903.8 


121.92 

400 

2.68 

8.8 

T-« 

245.09 

804.1 

G-9 

182.88 

600 

3.81 

12,5 

T-7 

214.95 

705.2 

G-10 

243-84 

800 

4.91 

16.1 

T-6 

184.04 

603.8 

6-11 

304.80 

1000 

6.16 

20.2 

T-5 

153.65 

504.1 

G-1 2 

365.76 

1200 

6.83 

22.4 

T-A 

123.06 

403.7 

G-13 

426.72 

1400 

7.86 

25.8 

T-3 

92.78 

304.4 

G14 

487.68 

1600 

8.14 

26.7 

T-2 

62.24 

204.2 

6-15^ 

-544.07 

-1785 

- 9.42 

-30.9 

T-1 

31.67 

103.9 

6-16** 

536.45 

1760 

7.16 

23.5 


‘Altitude of tower base is 1 1 1 2.5 m (3650 ft) MSL 

^hese microphones were positioned on a line normal to the flight path (see fig. 2). 


The BREN tower was also used as a platform for observing the structure of the lower 
atmosphere while the sonic boom measurements were being made. Anemometers and thermometers 
were located at eight levels. In addition, low-level measurements were made at two locations on a 
30-m (100-ft) mast located 200 m (656 ft) to the south of the BREN tower. This position was 
selected to avoid distortions caused by the small buildings nearby. The heights of the sensors are 
given in table 3. Figure 4 shows some of the temperature and wind sensor instrumentation attached 
to the tower by A-frame supports. The microphone supports can be seen extending from the right 
side of the tower. 
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TABLE 3.-L0CATI0N OF METEOROLOGICAL INSTRUMENTATION 
ON THE BREN TOWER 


Level no. 

Height of wind systems 

Height of thermometers 

m 

ft 

m 

ft 

9® 

460.3 

1510 

459.8 

1508 

8^ 

303.3 

995 

300.0 

984 

7 

217.9 

715 

214.6 

704 

6 

175.3 

575 

171.6 

563 

5 

132,6 

435 

129.0 

423 

4" 

89.9 

295 

86.3 

283 

3 

47.2 

155 

43.6 

143 

2 

25.9 

85 

24.4 

80 

1 

15.2 

50 

14.6 

48 

0 

3.1 

10 

2.4 

8 


^Levels so indicated contained a vertically-orientated propeller anemometer. 


A Piper Comanche airplane was used to obtain additional data on the atmospheric conditions 
above and in the vicinity of the BREN tower. Ljpper level temperature, wind, humidity, and 
pressure conditions to about 12.2 km (40 000 ft mean sea level) were obtained by using standard 
rawinsonde equipment at the Jackass Flats weather station, 4JA, located 3 miles northwest. 

Sonic booms were generated by NASA F-104 airplanes. The flight track was nominally on a 
true heading of 035® in line with the ground microphone array, except for the lateral cutoff flights 
which were on a 125° heading and displaced to the south by 18.4 to 24.2 km (1 1.4 to 15.0 st mi). 
The airplane altitude above mean sea level ranged from 9.78 to 10.30 km (32 100 to 33 800 ft). All 
flights were controlled by radar and the flight tracks were recorded for later analysis. Each test 
airplane made three to four consecutive passes over the BREN tower. Voice communication with 
the pilot provided the capability to change the flight conditions appropriately on the next pass 
depending on what was observed at the tower. 


12 



THEORETICAL SONIC BOOM CALCULATION INPUTS 


The primary objectives of the study of the BREN tower measurements included (1) an 
interpretation of the observed data guided by theoretical study of the test conditions and (2) a 
determination of the validity and accuracy of current sonic boom calculation methods. To 
accomplish these the computer method described in references 7 and 30 was used. Several types of 
data were required as inputs, which include the airplane flight conditions and the meteorological 
conditions between the airplane and the ground. Accordingly, the necessary data were recorded 
during the test series for later analysis. In addition, the airplane lift and shape characteristics were 
required in terms of the airplane F-function. This section contains a discussion of the methods used 
to obtain each type of input data. 


Flight Conditions 

Among the important flight conditions measured by tracking radar were the airplane altitude, 
location with respect to the BREN tower, and airplane velocity as a function of time. These 
variables were specified as accurately as possible, since relatively small inaccuracies can lead to 
unrealistic theoretical predictions, particularly for low supersonic Mach number flights. 

In most cases it was necessary to smooth the observed radar data, since the radar accuracy 
tolerance was significant. Figure 5 shows an extreme case of variation of radar-observed airplane 
ground speed with flight time for one flight. Tolerances in the radar accuracy and a ‘‘smoothed” 
curve are indicated on the figure for reference. The large change in ground speed between 23 and 28 
sec corresponds to an axial load factor (longitudinal acceleration) of about 1.7 m/sec^ (5.5 ft/sec^) 
or 0.17 g. This particular flight, however, was a “steady” flight near the threshold Mach number. 
Since most of the variation observed appears to be due to accuracy tolerance, the maneuver data for 
each case was smoothed for analysis purposes using a least-squares fit to the data. This smoothing 
assured reasonable analytical results and reduced the influence of recording accuracy. Similar 
techniques have been used by the French in analyzing comparable data (see the appendix of 
reference 26). 

In addition to the radar-observed flight conditions, pilot-read altitude, Mach number, and 
fuel-on-board were also available at both the steady and breakoff points. In general the pilot-read 
altitudes were used; the pilot-read Mach numbers were not used because of instrument inaccuracies 
at low supersonic Mach numbers (except in cases where no radar data were available). 


For most cases the airplane location and flight time were known, so that the theoretical boom 
arrival time at the tower could be calculated and compared with actual arrival times. This 
comparison provided a good test of the accuracy of linear geometric acoustics theory. 


Meteorological Conditions 

During the flights at low supersonic Mach numbers, the sonic boom traveled a relatively long 
distance through the atmosphere, so that meteorological conditions had an important effect. In 
general, sensitivity to atmospheric conditions increases with decreasing Mach number (refs. 31, 32 
and 33). Because of this, an effort was made in the BREN tower tests to obtain accurate 
measurements of the meteorological conditions, particularly in the lowest kilometer. Efefinitive 
meteorological measurements were made by using several data-gathering methods. These included at 
least two rawinsonde observations during the boom-generating periods, two to four Piper Comanche 
runs, and the BREN tower data at the time of each boom (see ref. 34 for further details). Thus, 
excellent measurements of the prevailing meteorological conditions in the test area were obtained. 

In the theory of references 7 and 30 a horizontally stratified atmosphere is assumed (variations 
in the vertical direction only). The effects of small-scale turbulence are ignored; only the relatively 
large-scale variations of temperature, pressure, and wind are taken into account. For each sonic 
boom overflight it was necessary to construct vertical profiles of temperature, pressure, and wind 
from the ground to the airplane altitude. This was done by using the BREN tower meteorological 
data (temperature and wind only), the Piper Comanche airborne turbulence measuring system 
(ATMS) temperature and pressure data, and the Jackass Flats weather station (4JA) rawinsonde 
observations of temperature, pressure, humidity, and wind. All of these data are given in reference 
34, along with a description of the instrumentation, recording, and analysis methods. Table 4 gives a 
sample data tabulation of a typical atmospheric sounding constructed from all sources. In compiling 
these data it was necessary to calculate pressure and humidity at the BREN tower levels based on 
the nearest rawinsonde observation in time, and to interpolate wind and humidity data at the ATMS 
significant levels, as noted. 

Table 4 also contains a parameter called “virtual” temperature, Ty, which is a corrected 
temperature used to account for the presence of water vapor. Appendix A contains a discussion of 
the calculation of the virtual temperature and its effect on shock propagation speed. Normally this 
effect is negligible, but for shock propagation near cutoff it may become important. Dewpoint 
temperatures from the rawinsonde observations were used as being representative of the conditions 
at the BREN tower. 
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TABLE 4.-TYPICAL METEOROLOGICAL DATA CONSTRUCTED FROM BREN TOWER, ATMS, 
AND RAWINSONDE DATA 
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‘Dewpoint temperatures and pressures were determined from rawinsonde data 
^Dewpoint temperatures and wind data were determined from rawinsonde data 


















































Significant meteorological variations can occur horizontally along the path that the boom 
travels that would not be measured by vertical soundings. This is particularly true for some of the 
threshold Mach number flights, when the boom passed within 1 .0 km (3300 ft) of the top of Little 
Skull Mountain (see ref. 34). 


Airplane F-Function Data 

A basic requirement for calculating theoretical sonic boom pressure signatures is the initial 
airplane disturbance, which is then propagated through the atmosphere to the ground. This was 
determined by calculating the F-function from the detailed geometry and lift distribution for the 
F-104 airplane. A brief description of the method used and tabulations of the resulting data at Mach 
1.3 and 1.1 are given in appendix B. The F-104 airplane is about 15 m (50 ft) long, has a reference 
wing area of 18.21 m^ (196 ft^) and, during the tests, had a gross weight between 7300 and 
9100 kg (16 000 and 20 0001b) depending on the fuel consumed. Four different F-104 airplanes 
from the NASA Edwards Flight Test Center were used. 
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I 

I 


DATA ANALYSIS METHODS 

f 

I 

i 

' In this section the methods used in the analysis of the meteorological and sonic boom data are 

briefly discussed. These include the calculation of the propagation speed from the meteorological 
' observations, the calculation of the shock propagation speed from the shock wave arrival times on 

the ground microphone array, the calculation of shock wave profiles from the tower shock wave 
arrival times, and the calculation of the ground reflection coefficient. 

I 

Acoustic Propagation Speed 

I 

i Directly beneath the airplane the propagation speed, Vp, is the speed of sound as determined 

from the temperature and wind speed in the direction of the normal to the shock wave. It is 
I calculated from the equation: 

Vp = a-u cos (i'-??) = a - Uj^ (1) 

The term u cos ( »> - t?) is the wind component in the direction of flight (tailwind negative), when the 
ray directly beneath the airplane is being considered. For the threshold Mach number and 
acceleration flights the airplane heading angle, was 035° (also v = 035°). For the lateral cutoff 
; flights, however, ^ was 1 25° while the shock wave heading angle, u, was about 90° at the BREN 

1 tower. The speed of sound is conveniently calculated from the temperature as follows: 

a = 20.045 [T(°C) + 273.2] j^/sec (2) 

In most cases where Vp is used in this report the sound speed, a, is corrected for the presence 
of water vapor by replacing the observed temperature, T, with a slightly higher temperature called 
I virtual temperature, Ty. This corrected sound speed is denoted by a* and the corresponding 

propagation speed by Vp*. A discussion of the calculation of the virtual temperature and its effect 
on the change in propagation speed is contained in appendix A. The propagation speed, Vp*, is of 
interest, since it is the speed at which an acoustic disturbance will propagate in the direction of 

flight. For the threshold Mach number flights and lateral cutoff flights V^* is of particular interest, 

\ ^ 

since cutoff will occur at the altitude where it equals the airplane ground speed, Vq. For the lateral 

cutoff flights it was necessary to account for the fact that the flight path was not parallel to the 

microphone array. 
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Shock Propagation Speed 


Since a two-dimensional array of microphones was used on the ground, it was possible to 
compute the velocity of propagation of the shock wave from the shock wave arrival times. For the 
threshold Mach number and acceleration flights where the flight track was parallel to the 
microphone array G-1 to G- 14, the shock propagation speed is: 

AX 

'^Pmic"lto(G-«4)-to(G-l)] 


where: 

AX = 975.4 m (3200 ft) 

) = shock wave arrival time at microphone ( ). 

For the lateral cutoff flights was calculated from the arrival times at microphones G-1 5 and 

, t^mic 

G-1 6 in an analogous manner. 

For cases where discrete shock waves were observed, . is a direct measure of the airplane 

t^mic 

ground speed, Vq. In many cases this value of the airplane ^ound speed was more accurate than 
the ground speed obtained from the tracking radar because of the accuracy tolerances discussed in 
the previous section. In some cases it was possible to note a change in this calculated propagation 
speed over the length of the ground microphone array. This was true not only for the longitudinal 
acceleration flights but also for some of the threshold Mach number flights. 

For cases where cutoff occurred above the tower and acoustic disturbances were observed 

(rather than shock waves), . is approximately equal to the propagation velocity, V^, 

Pmic ^ 

determined from the meteorological data. This was to be expected since Vp (or Vp*) is the 

propagation velocity of acoustic disturbances, which is dependent on the temperature and wind 

conditions near the ground. . was calculated for all cases where it was possible to identify 

pressure signature features from microphones G-1 to G-1 4. For some rumble cases it was not 

possible to make this calculation. 


Shock Wave Profiles 


The sonic boom measurements on the tower provided the capability to calculate the shock 
wave front shapes and locations at a reference time in the vertical plane. This required a conversion 
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from time to distance, since the shock wave front swept past the fixed tower at a given velocity, 

V- . . The shift in distance from the tower, AX, corresponding to differences in shock arrival 
.t'mic ^ . 

time, t^ from a reference time, tj^gp, is: 


AX-(tREp-to)- Vp^.^ (4) 

The reference time was taken as the arrival time at tower microphone T-1. Since the microphones 
were positioned vertically on the tower, no correction because of microphone offset was necessary. 
The analysis of shock wave geometries in the vertical plane added significantly to the interpretation 
of the observations and to comparisons with linear theory. 

An associated calculation is the conversion of the observed pressure signatures (overpressure 
versus time) from the time scale to a distance scale (overpressure versus distance). By making this 
conversion, shock wave profiles and observed pressure signatures could be placed on the same graph 
using a common distance scale. 


Ground Reflection Coefficient 

The ground reflection coefficient, Kr, was of particular interest for the BREN tower flight 
tests, since shock waves near the cutoff condition were observed. For shock waves near cutoff, 
incident on the ground, the reflection coefficient was expected to approach 1.0. (For oblique shock 
waves incident on smooth surfaces, the reflection coefficient is close to 2.0.) In this study Kr was 
determined by comparing intensities on the tower with those on the ground. The following 
equation was used: 


[ AP„,,,(G-7>tAP^,,(G-8)l 

" lAPmaxfr-l) + AP„^^(T.2)l <5) 

By using two microphones for both tower and ground maximum overpressures, the effects of 
small-scale atmospheric turbulence were minimized. The specific microphones selected roughly 
correspond to those that recorded the same portion of the shock wave during and after reflection 
from the ground. Even with this technique, however, some inaccuracy resulted since large 
overpressure variations can occur over small distances. In some cases it was necessary to examine 
pressure signatures in detail. 



ANALYSIS OF THE THRESHOLD MACH NUMBER 


FLIGHT TEST DATA 


This section contains the results of the analysis of the sonic boom characteristics measured 
during the threshold Mach number flights. The introduction contains background material on the 
nature of the phenomenon and the definition of the threshold Mach number. A summary table of 
the flight conditions for the 79 threshold Mach number flights is given. The presentation of the 
study results has been separated into categories according to what was observed on the BREN 
tower. 


Introduction 


Shock wave cutoff at low supersonic speeds .— criterion for shock wave cutoff above the 

ground from a supersonic airplane is that the airplane ground speed must be less than the maximum 

speed of propagation of the shock wave beneath the airplane. This concept is illustrated in figure 6. 

An airplane with ground speed, Vq flown close to the ground below altitude Zj will be flying 

subsonically, so that no shock wave will form. For flight at a higher altitude, Z 2 , however, with the 

same ground speed, the airplane will be operating supersonically due to the lower sound speed 

there, and a shock wave will form. In that case, the shock wave is '‘bent” or refracted by the shock 

propagation speed gradient as it travels toward the ground and will be completely refracted at 

altitude Zj where the propagation speed is equal to the airplane ground speed. At the altitude of 

complete shock wave refraction, the shock wave path is parallel to the ground. As the airplane 

ground speed or Mach number increases, the initial shock wave angle is less vertical with respect to 

the ground, so that it takes larger gradients of temperature and wind for the shock wave to become 

normal to the ground. For an airplane ground speed equal to in figure 6a, the shock wave 

t'max 

will be completely refracted at the ground (see fig. 6b). For an airplane ground speed greater than 
, the shock wave and a sonic boom will be observed at the ground (see fig. 6c). 

t’max 


At the altitude of complete shock wave refraction, some interesting phenomena occur. First, a 
slightly more intense shock wave may occur due to the fact that the energy is focused, producing a 
caustic. The initial shock wave intensity is relatively low, however, due to the fact it was produced 
at a Mach number slightly greater than unity. Below the altitude of complete refraction, the shock 
wave and sonic boom degenerate to acoustic or sound waves which propagate toward the ground. 
These acoustic waves are perceived as rumbles similar to distant thunder rather than the “bang” of a 
sonic boom. If the altitude of refraction is high enough above the ground, no noise will be 
discernible at the ground. 
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Threshold Mach number. —The threshold Mach number has been defined to be the maximum 
airplane Mach number for which complete shock wave refraction can occur at or above the ground. 
For cutoff directly beneath the airplane, the equation defining the threshold Mach number (ref. 32 
and appendix C) is: 


where: 


max } 


( 6 ) 


Mn 


= threshold Mach number 


= altitude 


a(Z) 


= speed of sound at altitude Z 


Un(Z) 


wind component at altitude Z parallel to flight path (tailwind 
is negative) 


[a(Z) - Uj^(Z)] = maximum shock propagation speed between the airplane 

and the ground in the direction of flight, with a(Z) and u(Z) taken 
at the same altitude 

Hq ~ sound speed at airplane 

Uj^ = wind speed at airplane (tailwind is negative) 


Equation (6) can be rewritten as: 


V<3=[a(Z)-u„(Z)|„,, = Vp^^^ (7) 

where: 

Vq = airplane groimd speed for flight at the threshold Mach number 

= (MT)(ao)-X 

[a(Z) - Ujj(Z)] = maximum shock propagation speed in the direction of flight 

Equation (7) simply states that for flight at the threshold Mach number, defined by equation 
(6), the airplane ground speed, Vq, is equal to the maximum shock propagation in the direction of 
flight speed, Vp^^^, beneath the airplane. More generally, the shock wave will be completely 
refracted at the altitude where the shock propagation speed is equal to the airplane ground speed. 


For steady, level flight near the threshold Mach number, the following statements can be 
made: 

• The ground speed of the airplane will be equal to the speed of sound at the cutoff 
altitude plus or minus the wind speed there, depending on whether the wind is a tailwind 
or headwind. 

• Although the allowable airplane speed is determined by the shock propagation speed at 
the cutoff altitude in the lower atmosphere, the allowable airplane Mach number is 
dependent on the temperature and wind at the airplane altitude. For a fixed ground 
speed, a headwind at the airplane gives an increased Mach number, while a tailwind results 
in a decreased Mach number compared to a case with no wind. 

Variability of threshold Mach number. -One of the potential advantages of commercial 
threshold Mach number operation is that the airplane ground speeds show little variation with 
direction of flight. Table 5 compares subsonic flight with threshold Mach number operation (where 
cutoff occurs at a safe altitude above the ground) over the San Francisco to New York City route. 
For average tailwind and headwind conditions the airplane ground speed varies by only about 
7.7 m/sec (15 kn), while for subsonic airplanes the variability is about 46.3 m/sec (90 kn). The low 
variability of ground speed with flight direction is due to the fact that the ground propagation 
speeds do not vary appreciably with season or flight direction. Table 6, compiled from the data in 
reference 35, illustrates average and extreme threshold Mach numbers, ground speeds, and block 
times over the San Francisco to New York City route. Further data on the variability of the 
threshold Mach number in time and space are given in reference 35. 

TABLE 5.-TYPICAL GROUND SPEEDS AND BLOCK TIMES FOR VARIOUS 
MACH NUMBER REGIMES OVER THE SAN FRANCISCO TO 
NEW YORK CITY ROUTE 

(Airplane altitude of 13.72 km (45 000 ft)) 



Mach 0.95 

Mach 0.83 

Threshold Mach number operation 
(Variable Mach number) 

Ground 

speed 

Block 

time, 

hr 

Ground 

speed 

Block 

time, 

hr 

Most probable 
route mean 
Mach number 

Most probable 
route mean 
ground speed 

Most 
probable 
block 
time, hr 

m/sec 

kn 

m/sec 

kn 

m/sec 

kn 

NO WIND 

280.4 

545 

4.6 

244.4 

475 

5.2 

1.15 

335.9 

653 

3.94 

TAILWIND, 
23.1 m/sec 
(45 kn)® 

303.5 

590 

4.3 

267.5 

520 

4.8 

1.08 

340.0 

661 

3.90 

HEADWIND, 
23.1 m/sec 
(45 kn)® 

257.2 

500 

5.0 

221.2 

430 

5.7 

1.20 

332.3 

646 

3.98 


® Mean wind component along the route at 13.72 km (45,000 ft) for the data used. 
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TABLE 6.-M0ST PROBABLE AND EXTREME ROUTE MEAN SAFE^ THRESHOLD MACH NUMBERS, 
GROUND SPEEDS, AND BLOCK TIMES FOR THE SAN FRANCISCO TO NEW YORK 
CITY ROUTE 
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'A safety factor consisting of a 5.1-m/sec |10-kn) reduction in the maximum airpiane ground speed 
was used, since cutoff must occur above the ground to avoid objectionabie sonic boom noise. 


Previous flight tests,— There; have been two experimental tests during which threshold Mach 
number flights were conducted. The first experiment (ref. 14) involved seven flights all on one day; 
the second (ref. 18) involved five flights, also on one day. Both experiments indicated that for 
airplane Mach numbers less than the threshold Mach number, neither boom nor rumbles were 
observed at the ground, whereas, for flight at Mach numbers greater than the threshold value, 
definite shock waves were observed in the pressure signatures. These earlier experiments, however, 
lacked the required density of microphones and number of flights required to define the shock wave 
characteristics associated with complete refraction or the cutoff condition. 

Effect of Terrain Slope 


One of the criteria for determining nearness to complete refraction was the angle of the shock 
wave front with respect to the ground. Complete refraction occurs when it is normal to the 
horizontal ground. This angle could be computed from the shock arrival times on the tower. During 
the analysis of the shock shapes it was noted that the angle of the incident shock wave with respect 
to the tower was not equal to the angle of the ground-reflected shock wave. Typically, the reflected 
shock was at a larger angle. This effect is attributed to the sloping terrain uptrack of the tower. 
Figure 7 illustrates schematically the effect of the sloping terrain. The reflected shock wave is 
rotated backward 2 «q degrees because of the degree ground slope. Two flights were chosen to 
verify the sloping terrain effect. These flights were at high enough Mach numbers to minimize the 
effects of changing airplane ground speed and atmospheric variations. In each case the deviation of 
the reflected shock front from being symmetric with the incident shock was determined with 
altitude on the tower. These data are given in table 7. 

TABLE 7.-CALCULATI0N OF EFFECT OF GROUND 
SLOPE ON REFLECTED SHOCK FRONT 


Tower 

mic 

no. 

Altitude 
above 
tower 
base, Zq 

Pass 026; 8-25; bongo 4-2 

Pass 082; 10-23; bongo 1-1 

Deviation of 
reflected 
shock from 
being sym- 
metric, AX 

Calculated 
ground slope, 

«G = 

1/2 tan-1 (^x/Z) 

Deviation of 
reflected 
shock from, 
being sym- 
metric, A X 

Calculated 
ground slope, 

1/2 tan'l ( AX/Z) 


m 

ft 

m 

ft 

rad 

deg 

m 

ft 

rad 

deg 

15 

458.4 

1504 

18.7 

61.4 

0.02036 

1‘10' 

19.6 

64.3 

0.02123 

1*13' 

14 

428.1 

1405 

17.9 

58.7 

0.02094 

1® 12' 

16.8 

55.1 

0.01978 

f 8' 

13 

397.4 

1304 

16.4 

53.8 

0.02065 

i“ii' 

17.2 

56.4 

0.02182 

1*14' 

12 

366.9 

1204 

14.9 

48.9 

0.02036 

f 10' 

14.3 

46.9 

0.01949 

f 7' 

11 

336.5 

1104 

13.7 

44.9 

0.02036 

fio' 

13.6 

44.6 

0.02036 

f 10' 

10 

306.0 

1004 

12.1 

39.7 

0.01978 

f 8' 

12.3 

40.4 

0.02007 

f 9' 

9 

275.5 

904 

11.7 

38.4 

0.02123 

1®13' 

10.5 

34.4 

0.01891 

1® 5' 

8 

245.1 

804 

9.7 

31.8 

0.01978 

1*8' 

9.0 

29.5 

0.01832 

1®3' 

7 

215.0 

705 

7.8 

25.6 

0.01803 

1® 2' 

6.6 

21.7 

0.01542 

0*53' 

6 

184.0 

604 

6.2 

20.3 

0.01687 

0*58' 

5.9 

19.4 

0.01600 

0*55' 

5 

153.7 

504 

- 

- 

- 

- 

4.6 

15.1 

0.01513 

0*52' 

4 

123.1 

404 

- 

- 

- 


4.0 

13.1 

0.01629 

0*56' 

3 

92.8 

304 


— 

“ 

— 

3.5 

11.5 

0.01891 

1*5' 
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The terrain slope with respect to the tower base between ground microphones G-1 and G-7 
uptrack of the tower is about 0.02036 rad (1® 10’). Thus, the ground slope angle, ocq, calculated 
from the nonsymmetry of the reflected shock, agrees very closely with the actual ground slope to 
within a small fraction of a degree. As a result of this analysis, the incident shock wave angle was 
used to determine nearness to cutoff, because of the influence of the sloping terrain on the reflected 
shock wave angle. 

Summary of Flight Conditions 

The 79 threshold Mach number flights have been summarized in table 8 for reference and a 
unique flight number was assigned to each pass for identification. The airplane altitude and Mach 
number are values read by the pilot at the steady point, and the airplane gross weight was 
determined from the steady point fuel-on-board. The maximum observed overpressure of all the 
ground and tower microphones is given, along with subjective observer comments of the boom 
character. Each pass was also categorized according to the subjective observations of the boom 
character. These categories are listed below: 


THRESHOLD MACH NUMBER CATEGORIES 


Category 

Reported Subjective Observation 

Flight Condition 

R1 

No Rumble 




R2 

Low Rumble 

M 

< 

Mj 

R3 

Moderate to Heavy Rumble 




R4 

Light Boom 

M 


M-p 

B1 

Near Cutoff (0 < 0.01 745 rad ( 1 .0 deg)) 

M 


Mj 

B2 

Normal Boom (0> 0.01745 rad (1.0 deg)) 

M 

> 

My 

B3 

Acceleration Effects 

M 

> 

M-p 


Mach Number Greater Than Threshold Value 

Thirty-one threshold Mach number flights were at Mach numbers sufficiently greater than the 
threshold Mach number to produce well-defined sonic booms on the ground. Hence, the airplane 
ground speed was greater than the maximum shock propagation speed, Vp* , between the 

^ indite 

airplane and the ground. For these cases the airplane ground speed exceeded Vp* by a range 
from about 1 to 36 m/sec (3.3 to 1 1 8 ft/sec). 
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TABLE 8.-SUMMARY OF THRESHOLD MACH NUMBER FLIGHT TEST CONDITIONS 


Category 

f- CM 

a: oQ Qc 

CO CM CM 
OQ CC CO 

cc 

CM CM CO r-** 
QC CC QC CQ 

QC 

CO CM 
CO GO CQ 

CO CO CM 
cn DQ OQ 

CO CM CM CO 
OQ OQ QC OQ 

CM CM ^ 
OQ QC QC 

CM CM CO CM 
QC OQ QC OQ 

CM CM CM ^ 
CQ CQ DC 00 

0 

® t 

> c 

1 

CO c 
c 

) 

0 

1 

3 

5 

3 

3 


No Boom or Rumble 
Boom, Near Caustic 
Very Light Rumble 

Triple Boom 
Very Low Rumble 
Boom, Rumble 

Long, Low Rumble 
Long, Low Rumble 
Moderate Rumble 
Boom, Rumble 

Rumble, Boom 
Double Boom 
Boom, Rumble 

Rumble 
Double Boom 
Double Boom 

Near Caustic 
Double Boom, Rumble 
Light Rumble 
Boom, Rumble 

Rumble, Boom 
Light Rumble 
Moderate Rumble 
Double Boom, Rumble 

Rumble 

Double Boom, High Peak 
Mod-Heavy Rumble 
Triple Boom 

Double Boom 
Double Boom 
Rumble 

Double Boom, Rumble 

03 

1— 

3 

03 

Mic No. 

IX) 

1 

op CO 

6h6 

lO LO o ^ 

T— *— <— 

HHHH 

CM ^ ^ 
6 66 

G-15 

G-7 

G-5 

T-14 

T-8 

G-5 

G-14 

in 

q CO 
6 6 H H 

in in 'd- 

00 r- ^ 

6 6 H 6 

■d CD 

q 00 t- 

6666 

Maximum 
ved overpr 

CM 

H— 

1.060 

0.129 

0.941 

0.094 

1.42 

CM o 

o 0 ) in CM 
o CM 

CD CD d CD 

0.744 

1.49 

0.644 

0.381 

1.40 

0.956 

1.72 

1.295 

0.191 

1.20 

1.40 

0.146 

0.544 

1.32 

CD 

CD CM 00 P^ 
CM ^ ^ CD 
d CM d ^ 

1.23 

0.886 

0.186 

1.05 

v> 

O 

CM 

2 

50.75 

6.18 

45.05 

4.50 

67.99 

00 o ^ r>. 

00 If) CO CO 
^ CM ^ 

^ CO 

35,62 

71.34 

30.83 

18.24 

67.03 

45.77 

in o in CD 
CO q 1 -; ^ 
CM CM O) P^ 
C30 CD ID 

CO CD in o 

O CD O CM 
P^ CD CD CO 
CD CM CD 

12.74 

115.87 

23.17 
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A direct measure of nearness to cutoff is the shock wave angle of incidence with respect to the 
vertical,©. This can be calculated from the arrival times of the shock wave at the various tower 
microphones after a transformation is made from time increments, AT, to distance increments, AX. 
This conversion method was previously described in the “Data Analysis Methods” section. For the 
cases discussed in this section, ©ranged from 0.01745 to 0.419 rad (1® to 24®). Shock wave profiles 
and pressure signatures measured on the tower are presented and discussed. 

Comparison of theory with observation. —Table 9 summarizes the pertinent observed data for 
these cases, including the airplane ground speed, Vq the shock propagation speed, Vp*, the shock 
wave inclination angle,©, the average observed maximum overpressure on the tower, AP t„„„, and 
the shock wave arrival time, t^, at microphone T-1 . For these cases it is of particular interest to 
compare the observed data with theoretical results. Theoretical values of ©, AP, and t^ are also 
given in table 9. 

The error of the theoretical arrival time is a good measure of the accuracy of the theoretical 
calculations. For the August flights this error is about ±1.0 sec. For the October flights, however, 
much larger errors were calculated (about ±5.0 sec). The significantly larger errors during the 
October flights are due to the meteorological conditions. The gradient of propagation speed in the 
lower 1 km (3300 ft) was nearly zero or negative during most of the flights. For this condition, 
small inaccuracies in the flight or atmospheric conditions would lead to large variations in 
theoretical shock locations and arrival times. Determination of the airplane ground speed is 
especially sensitive. For example, starting with a ground speed 3.5 m/sec (1 1.5 ft/sec) greater than 
Vn* , a 1 .5 m/sec (4.9 ft/sec) decrease in the airplane ground speed will change © from 8® to 6°, 
with accompanying changes in AP and t^. An additional 2 m/sec (6.6 ft/sec) decrease in ground 
speed will change 0 from 6® to 0®. Thus, the accuracy of the theoretical calculations is highly 
dependent on the accuracy of the measured airplane ground speed, as well as the accuracy of the 
available meteorological data in representing the conditions along the boom propagation path. 

Figure 8 summarizes the comparison between theory and experiment for the shock inclination 
angle,©. The values of ‘ Vp*) ^ table 9 and the sound speed, a(Z), near the ground were 

used to calculate theoretical values of ©from appendix equation (C8), which is: 



A convenient relationship between (Vq . - and (M - M j) is given by appendix equation 

(CIO) as: 


(M-MT) = (VG-Vp*^J/a„ 
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TABLE 9.-SUMMARY OF THRESHOLD MACH NUMBER CASES FOR WHICH M > M j (CATEGORIES B2 AND B3) 
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'Calculated from shock wave arrival times at ground microphones 
^Average value between 100 and 500 ft on tower 
'Equation (C8I was used with (Vq - Vp*) values 



























TABLE 9-SUMMARY OF THRESHOLD MACH NUMBER CASES FOR WHICH M > Mj (CATEGORIES B2 AND B3) 
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Thus, an (M - Mt) scale is also given in figure 8. The observed (Vp > V„* 

^ ^ max 

reasonably well with theory. 

Examination of the overpressure data in table 9 shows that in many cases relatively high values 
occurred for flight at Mach numbers well above the threshold value. These cases typically exhibited 
caustic-like pressure signatures, which were apparently produced by relatively low-magnitude 
airplane accelerations rather than by ray focusing due to the cutoff condition. These cases have 
been placed in a separate category, B3. Relatively low-magnitude accelerations from below the 
threshold Mach number could produce caustics similar to accelerations begun at subsonic Mach 
numbers. Variations of 5 m/sec (16.4 ft/sec) or more in the airplane ground speed were noted from 
the radar data, and these would be sufficient to produce caustics. In some cases the theoretical 
calculations are influenced by changes in the radar-observed airplane ground speed. In general, an 
acceleration results in overpressures higher than steady flight, while a deceleration results in lower 
overpressures. During steady flight at about Mach 1.2 the nominal theoretical maximum 
overpressure is about 28.7 N/m^ (0.60 Ib/ft^) for shock waves observed on the tower (this value is 
approximately doubled at the ground due to reflection). 

Figure 9 shows the variation of maximum overpressure measured on the tower near the cutoff 
condition. The cases where caustics were produced by accelerations are indicated by triangular data 
points. The cases for which cutoff occurred on the tower are plotted at zero (Vq - Vp*) and are 
discussed in detail in a later section. These data indicate that the linear theory is fairly accurate to 
within (M-Mj) = 0.01 and (Vq - Vp*) = 2.0 m/sec (6.6 ft/sec). In addition, these data indicate 
that accelerations can have an important effect near cutoff. 

In general, the agreement between the theoretical calculations and the experimental data given 
in table 9 is good. The range of differences between the theoretical calculations and the 
observations for©, AP, and t^is as follows: 

e «« ±0.0436 rad (2.5°) 

AP (with respect to APj) ±17.0 N/m^ (0.35 Ib/ft^) 

t^ ±1 .0 sec (except ±5.0 sec for cases with inversions) 

In view of the sensitivity of the theoretical calculations to the airplane ground speed, the airplane 
ground speed as determined from the shock propagation speed over the ground microphones, 
VQmic? was used instead of the ground speed measured by the tracking radar for several cases. 
These results are given in table 10 for passes 015, 054, and 109. For these three cases use of Vq^.^ 
gives substantially better agreement. 
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TABLE 10.-EFFECT ON THEORETICAL BOOM ARRIVAL TIME OF 


USING Vq . INSTEAD OF RADAR GROUND SPEED 
'^mic 


Date 

Bongo- 

pass 

Pass 


error using 
sec 

t error from 
table 9, 
sec 

m/sec 

ft/sec 

8-25 

1-2 

015 

364.1 

1194.5 

0.0 

+0.4 

8-31 

1-4 

054 

354.2 

1162.1 

-0.1 

-0.5 

10-29 

2-2 

109 

344.6 

1130.6 

I 

1 

-2.4 


Shock wave profiles and pressure signatures.— Tht shock wave profiles and tower pressure 
signatures are shown in figures 10 through 13 for passes 101, 106, 1 17, and 028, respectively. The 
incident and reflected shock waves (onset of each) are plotted to scale so that © values can be read 
directly. In addition, the pressure signatures were converted to the same distance scale and are given 
along with the shock wave profiles. A time scale and overpressure scale are also shown for the 
pressure signatures. 

The data for pass 101, given in figure 10, is a case where acceleration effects were not of great 
importance. The angle 0 for this case is about 0.105 rad (6“). A theoretical pressure signature is 
superimposed on the observed pressure signature at microphone T-13; the agreement in this case is 
very good. The incident N-wave at the lower tower levels appears to have been distorted by 
small-scale atmospheric turbulence. 

As noted in the previous section, in several cases accelerations produced caustics for cases 
where the incident shock wave was not near cutoff. Two such cases are shown in figures 1 1 and 1 2 
(passes 106 and 117, respectively). In both these cases a decayed pressure wave or pressure pulse 
precedes the basic caustic-like pressure signature. The angle of the incident shock front with the 
vertical,©, is about 0.1745 rad (10^. 

Pass 028, shown in figure 13, is a particularly interesting case since © at the ground is about 
0.01745 rad (l.(f). Here again, however, it appears that acceleration effects have produced the 
caustic-like pressure signatures rather than focusing due to the cutoff condition alone. 

In several cases a series of pressure signatures were observed. In one case (pass 004) at the 
ground a “triple boom” was observed. Changes in airplane ground speed (decelerations as well as 



accelerations) can produce such effects, since the propagation path is very dependent on the 
airplane Mach number, and shock waves emitted from different portions of the flight path may 
reach the same point at the ground at about the same time (see fig. 5). 


Mach Number Near Threshold Value 

Eleven flights produced shock waves that were cut off on the tower. The shock wave 
inclination angle for these cases was between 0 and 0.01745 rad (0 and T). The airplane ground 
speed was thus within 0.1 m/sec (0.3 ft/sec) of the shock propagation speed on the tower, and the 
airplane Mach number was within 0.0005 of the threshold Mach number. For these cases it was not 
possible to perform theoretical calculations, since the shock waves were too near cutoff. Table 1 1 
contains a summary of the observed data for these cases. 

It is interesting that the overpressure intensities on the tower for the cutoff cases are not 
appreciably greater than would be produced during steady flight at higher Mach numbers (about 
28.7 N/m^; 0.60 Ib/ft^), with the exception of pass 063 (see table 1 1). Thus, it would appear that 
accelerations through the threshold Mach number can produce caustics stronger than those 
produced by the cutoff phenomena. 

Precursors. —k characteristic of pressure signatures near the cutoff condition is that pressure 
pulses or “precursors” are frequently evident propagating ahead of well-defined shock waves in the 
pressure signature. Precursors appear to be associated with near-sonic conditions, where disturb- 
ances generated by the shock wave can propagate ahead of it. In some cases several precursors could 
be seen. Henderson (ref. 36) has presented some theoretical work on this aspect of the cutoff 
phenomena. He shows that precursors can be stationary with respect to the shock wave (“bound”) 
or nonstationary (“free”). Bound precursors may be thought of as a detached shock moving ahead 
of a blunt body. The author shows that free precursors may be produced before actual cutoff. 
Postcursor waves propagating in the downstream direction also occur, but these were difficult to 
differentiate in the experimental data from acoustic-like disturbances that normally trail behind the 
main part of the pressure signature. In Figures 14 through 20 the onset of precursors is indicated by 
dashed lines and the onset of shock waves by solid lines. 

Shock wave profiles and pressure signatures. — 14 through 20 contain shock wave 
profiles and pressure signatures for the most interesting cases, where cutoff occurred on the tower. 
The pressure signatures are given to the same distance scale as the shock wave profiles; a time scale 
is also given for the pressure signatures. 
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TABLE 11.-SUMMARY OF CASES FOR WHICH CUTOFF OCCURRED NEAR THE GROUND (CATEGORY B1) 
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Passes 002 and 114, shown in figures 14 and 15, respectively, are two cases where cutoff 
apparently occurred very near the top of the tower. In both cases caustic-like pressure signatures 
occurred at the tower top preceded by a pressure rise, coincident with the shock wave becoming 
vertical to the ground (0=0). The shock waves decay rapidly with distance below the cutoff 
altitude, with acoustic-like disturbances occurring at the ground (but still identifiable as decayed 
pressure signatures). 

For three cases cutoff occurred near midtower. These are shown in figures 16, 17, and 18. In 
each case caustic-like pressure signatures occur simultaneously with the shock wave becoming 
vertical (0=0). Pressure pulses propagating ahead of the basic signatures are evident. Passes 063 and 
102 (figs. 17 and 18) are particularly interesting, since the incident and refracted shock waves that 
form the caustic at the cutoff altitude are evident. The pressure signature observed at microphone 
T-11 during pass 063 had the largest overpressure measured during the threshold Mach number 
flights. This maximum value was 135.0 N/m^ (2.82 Ib/ft^). This pressure signature, however, 
appears to have been influenced by microscale turbulence. 

Pass 103 (fig. 19) shows the effect of the varying propagation speed on the shock front 
orientation. In this case cutoff may have occurred above the tower because of an inversion in the 
propagation speed gradient, since 0 is negative at the tower top. Caustic-like pressure signatures are 
evident near the tower base where a second complete refraction occurred. 

Overpressure variation near caustics. -The variation of shock wave strength in the vicinity of 
the caustic was determined from the data. Figure 21 gives the observed maximum overpressures 
with distance from the caustic for several cases where caustics associated with the cutoff condition 
(0= 0) were observed on the tower. Amplified intensities occur only within about 100 to 200 m 
(330 to 660 ft) above and below the caustic. This suggests that linear theory is capable of predicting 
intensities to about 100 m (330 ft) above the caustic with reasonable accuracy, since the nonlinear 
effects that predominate at caustics become important at that distance. 

The maximum overpressure observed on the tower was 135.0 N/m^ (2.82 Ib/ft^) during pass 
063. This signature, however, appears to have been influenced by microscale turbulence, with the 
high overpressure due to “spikes” superimposed on the front and rear shocks (see fig. 17). The next 
highest overpressure on the tower for these cases was about 51 N/m^ (1.06 Ib/ft^), and is more 
representative of the caustic intensity produced by threshold Mach number flight. Compared to 
steady, level flight at about Mach 1 .2 (which gives an intensity of about 28.7 N/m^ (0.60 Ib/ft ) on 
the tower), the caustic amplification is about 1.8. 

A fundamental consideration is whether the caustic intensity is related to the gradient of the 
shock propagation speed. A stronger caustic intensity would be expected where a large gradient in 
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propagation speed existed, and focusing would occur over a small altitude range. In this case there 
would be a rapid decay of intensity below the caustic. Passes 063 and 102 (figs. 17 and 18) can be 
considered to represent this situation. When the propagation speed gradient is small, however, the 
focusing would occur over a larger depth of the atmosphere, and only moderate amplification 
would be expected. In addition, caustic-like pressure disturbances may propagate a considerable 
distance below the actual caustic, particularly if an inversion exists in the propagation speed profile 
below the cutoff altitude. Figure 20 shows the data for pass 112, where cutoff occurred above the 
tower but caustic-like signatures of moderate intensity are observed on the tower. Thus, the 
meteorological conditions which produce the cutoff shock wave and associated caustic have a 
fundamental influence on the nature of the observed pressure signatures in the vicinity of the 
caustic. 


Mach Number Less Than Threshold Value 

Of the 79 threshold Mach number flights, 37 produced shock waves which were cut off above 
the tower. Relatively low-intensity acoustic disturbances were recorded on the tower and the 
ground, with overpressure levels less than 24.0 N/m^ (0.5 Ib/ft^). In general, observers subjectively 
rated the disturbances for these 37 flights as rumbles (low, moderate, or heavy). In a few cases, 
however, they were rated as “low” or “very slight” booms, and in five cases cutoff occurred at an 
altitude high enough above the ground so that no booms or rumbles were observed. These 37 
“rumble cases” provide information on the predictability of the cutoff altitude and on the nature of 
the acoustic noise that occurs below the cutoff altitude in terms of the known meteorological 
conditions. To aid in this analysis a number of measured and calculated parameters were determined 
and compared. These are shown in table 12. 

The first of these parameters is the calculated propagation speed of the disturbance over the 
ground microphone array determined from the arrival times at the microphone locations. An 
average value over the length of the ground array was used. This propagation velocity is called 
^Pmic* second parameter is the propagation speed, determined from the observed 

meteorological conditions in the vicinity of the tower. The effect of the presence of water vapor 
was taken into account (see appendix A). The maximum value of Vp* from the tower 
meteorological measurements was used. In some cases, however, a significant inversion occurred in 
the Vp profile (on October 23 and October 30) so that the maximum value above the tower was 
noted in the table instead, since it determines the cutoff condition. The airplane ground speed, Vq, 
was obtained from the radar data. Since the exact portion of the flight path which produced the 
cutoff shock wave above the tower was not known, the accuracy of this speed is about ±2 m/sec 
(±6.6 ft/sec). Other parameters in table 12 include the theoretical cutoff altitude, Z^, the average 
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TABLE 12-SUMMARY OF RUMBLE CASES-THRESHOLD MACH NUMBER FLIGHT 


Ground 
reflection 
coefficient, Kp 

1 00 00 LT> O 00 
' o> 00 ro o> CM 
d o — d 

00 r»v m CO 

CM 00 00 ^ o 

00 O) ^ CD 

CD p p ^ ^ 

1.66 

1.41 

r>' CM 
CO O CX) 

d — d 

CM CO r«. . . 

CD CD C- 1 1 

d d d 

^ 00 ^ ID 
r* Tt 00 00 

d d 

^ ID CO CM 
p p p P 

Cate- 

gory^ 

CN CM <N (N CO 

OC cr CC QC QC cc 

CO CM CM rr CM CO 
OC QC QC QC OC OC 

CM ^ ^ ^ 

OC GC CE QC QC 

^ CO 

QC QC QC QC 

CO CM Tf 

QC QC QC 

CM CO CM ^ ^ 
QC QC QC QC QC 

Tf CM ^ CM 

QC QC QC QC 

^ CM ^ ^ 
QC QC QC o: 

xeu 

Ib/sq ft 

O) 00 — CM CO 

lO CO CD CN 

o o o o o »- 

d d d d d d 

) 

CM C3> g? ID o 

o ’cr o CM o tt 

CM T-’- CM CM CO 

d d d d d d 

(O 00 00 00 
^ CM CD ID O 
CM CO CO 

d d d d d 

CD CO 

CD 

0-0 0 — 
d d d d 
1 1 

0.175 

0.283 

0.403 

CO CO 05 
CD CM 

CM t- o O 

d d d d d 

CO CD O CD 
00 »- CD CO 
CM CM »- «- 

d d d d 

T- oo CO CM 
CO ID CD P 
CM T- ^ CM 

d d d d 

Is 

c\ 

E 

z 

CM CD O CO 

O 00 CM ^ O) p 

d CM CO CO (N CD 

1 

f^CO<N*-00»- 
<0 CM p p P 

o> ID d o> d 

Tf to CD CM 00 
CD CO ^ ^ 

CD d 00 ^ 

CO r- 

r- o O ^ 
a> d o ai 
^ 1 1 

00 CD 00 
to CD p 
00 CO 00 

p p p p p 

d 00 d d 
^ ? 1 

CO 05 CM CD 
CO ^ ID CD 

00 d to 

o o ^ 

ID p p p 
CM CO 

t- CM *- 

Q> 
0> "O 

« 

o o o o o o 
o o o o o o 

T- 00 ID ID m 

CD ^ ID CM CM »- 

888888 
CD o CO o 00 

CM CM CM CM CM 

o o o o o 
o o o o o 

O CO O CD CO 
CM CM CM •- CM 

CM. » 

888 

CO O CD 
— CM — 

o o o o o 
o o o o o 

05 CM CD CM CM 
ID CD ID ID 

o o o o 
o o o o 

ID CD ID 
ID CM 

8888 
CM CM CM CM 
CD rv 

m ^ rvj 

E 

o o o o o o 

CD CD CD CD CD CD 

O O O O ID O 
CD O O O CM ID 
^ CO CO 00 00 

o o o o o 
o o o o o 


O O O 
O O CD 

O O O O O 
O O O O O 
ID CD O CD CD 
CM •- r- 

o o o o 

ID p O p 

8888 



CD ^ CD ^ 



^ CO ^ 


•- CM CM CM 

u 

tj Q if 

« 

O O O O O O 

O O O O O O 

CO CO CO CO CO CO 

CD CO ^ ^ 

o o o o o 

O O O O O 1 
— CO 0> CO CO ' 
^ CM CO 

o o o o o 

o o o o 

<D O CD 

CO ^ «- 

1 O O 1 
CM CO 

7T 

o o o 
o o o 
00 
CM 

O O O Q O 

O o o o o 

CO CD O O CM 
CD ID 05 O 05 

o o o o 
o o o o 
P O 00 00 
^ CM CM ^ 

4900 

-8200 

^ 3 2 
c o .t; 
Q- H 

m 

E 

CD O CD O O O 
CM O CM O O O 
CM O CM O CO CD 
CM CM CM t- »- 

o o o o o 

CD O 00 O O 1 
CO CO 00 ^ o ' 

o o o o o 

CD 00 O CD 

CD CD CO ^ 

o o 
1 S O 1 

CM CM 

o o o 

CO to CD 
CO CM 

O O ID O O 
O ^ ID O 
O o CO 

CM »- CM CO CM 

o o o o 

O O ID ID 
ID CD 00 ID 

o o 

ID ID 1 1 

«- CM 
) 

lb 

> 

X 

— 

ft/sec 

CO CO CM CD O p 
d CO 00 

CM CM CO 

00 p r** p CO 

CO ^ ^ 

^ CM CM ^ CM , 

P CD p p p 

d <j> r^.' d d 

CO CM fM. CO 
CD CD CO 
CM «- 

^ p p 
cb 

»- CM 

^ ^ ^ CO 

^ CM CO 

t— CO Tt ^ 

05 CD — p 
^ CD CD 
CM t- CM 

^ 05 'O' 
CO lb CD lb 

t- •“ 

X3 T3 -O 

E 

* a 
> 

m/sec 

^ CM CD CD CD 

K CO CT) CD CO 

P CM P P KT 
<r r*^ CO CD d 

CO CD ^ O CM 
CO CM CM d CO 

'T p P p 

d 00 ^ ^ 

P 53- ^ 

cb CM 

CO 53- p o P 
lb 53^ d CO CO 

^ p p 

CM tx> 

p p p 

-«cr ^ CM 
T) "O T3 -D 

0) ^ 
c “P > 
«3 C ^ 

ft/sec 

1119.1 
1131.9 

1120.1 

1143.0 

1140.1 
1145.0 

1149.9 
1136.1 

1143.0 

1151.9 

1147.0 

1169.9 

1145.0 

1148.3 
1149.9 

1160.1 

1161.4 

p p p p 
CM CM CD 

^ CM to ^ 

1066.9 

1055.1 

1076.1 

^ ^ p ^ p 
^ d CO lb d 

CD ID CD 

o o o o o 

1100.1 

1113.8 

1113.8 

1112.8 

1129.9 

1135.2 

1145.0 

1145.0 

w O 0> 

■< 0)g_ 

m/sec 

<f~ O ^ p o 
CD 00 r*^ d 

^ ^ ^ ^ ^ 

CO CO CO CO CO CO 

350.5 

346.3 

348.4 
351.1 

349.6 

356.6 

p p p p p 
d d d CO ^ 

^ CD ID CD CD 
CO CO CO CO CO 

p p p p 
d CM CD d 
O' -O’ o- o- 
CO CO CO CO 

325.2 

321.6 

328.0 

^PppT- 
00 ^ CO 

CM CM CM CM CM 
CO CO CO CO CO 

p p p p 

lb d d d 

CO CO CO CO 
CO CO CO CO 

Tt o p p 
^ CD d d 
^ ^ ^ 
CO CO CO CO 

jrved 

near 

ax 

und 


CM CM p ^ p 
CO CM «— CD 00 CO 

00 p p p 

^ ^ Oi *— 00 

P p P P 

CD d »- 

CO 

c6 

1082.0 

1082.7 

1084.0 

p p TT P P 
^ CO CO 

p ^ p ^ 

r«i 00 d d 

p o 

d CO CO d 

£ 

^ CD CD CD 

(o CD CO CO CO 



05 05 05 05 O 

o o o o *- 


u u u 

irt c o 

S-. - 

u 

o 

CD CM C7) CD O CD 

o6 ^ d c\i CO CM 

^ CD CD CD CD CD 
CO CO CO CO CO CO 

o p p p p p 

CD CO cb CO rv.’ CD 
CD CD CD CD CD CD 
CO CO CO CO CO CO 

P P p P CM 
CM CM CM to 
CD CD ID CD CD 
CO CO CO CO CO 

p 

329.8 

330.0 

330.4 

CO ID CD CO 
CM CO ^ CD 

p p ^ p 

p p p p 

lb 00 00 d 

Tt ^ TJ- 

CO CO CO CO 
O U L> 

> 

E 

ID 

CO 

CO CO CO CO CO 
CO CO CO CO CO 

CO 5 

CO CO CO CO 

i L- o 

ft/sec 

No Data 

1 1 59.4 
1154.8 
1154.8 
1158.1 

1 1 58.5 

1163.7 

(b) 

1168.0 

1168.3 

1165.3 
1158.5 

1158.1 

1161.1 

1164.0 

1162.1 
1172.9 

1145.7 

(b) 

(b) 

1147.3 

1074.5 

1081.4 

1079.7 

(b) 

1100.1 
1102.0 
No Data 
No Data 

1102.3 
1136.1 
1128.9 

1143.4 

1143.0 

1150.6 

1151.9 

1160.4 

— O) a> F 

J 

ro O ^ 

m/sec 

<o 

Q^pOO»- 
_ CO CM CM CO CO 
O CD CD CD CD CD 
Z to CO CO CO CO 

354.7 

(b) 

356.0 

356.1 

355.2 
353.1 

p p p p p 
CO CO ^ ^ 

CD CD CO ID CD 
CO CO CO CO CO 

349.2 

(b) 

(b) 

349.7 

327.5 

329.6 
329.1 

(b) 

335.3 
335.9 
No Data 
No Data 

p p ^ p 
CD CD OC) 

CO ^ 

CO CO CO CO 

348.4 

350.7 
351.1 

353.7 

Pass 


CO CD r«. 00 CD 

O O o o o d 

O O O O O O 

Tf CD CM CO CD 
•- »- CM CM CM CM 
O O O O O O 

CO ID CD CD Tf 
CD ID CD cD CO 
O O O O O 

CM CO ^ CD 

00 m 00 00 
o o o o 

089 

090 

091 

CD r-* t» 05 o 

05 05 05 05 O 

O O O O <■- 

107 

111 

112 

115 

00 O »- CM 
«- CM CM CM 

Bongo- 

pass 

CO CM CM CO 

^ CM CO CO CO 

r- p p p r- p 

CM CO CO tT ^ 

p r- CM p 

CM CM CO 

•- CM CO ^ 

^ p p 
CM CM CM 

p CO p p 

CM CM CM 

^ tT ID CO 
CM CM CO 

CO «- CM CO 
CM CM CM 

Date 

8-24 

8-25 

8-31 

10-23 

10-27 

10-28 

10 29 

10-30 


E 

o 

o 

XJ 


E 


- o 


-O 

o 

E 




o 


o ^ 


^ ^ 0 ) 


Q 

x: 


38 



observed maximum overpressure of all 16 ground microphones, subjective category of 

ITlaX 

the disturbance, and the ground reflection coefficient, as determined from several of the ground 
and tower maximum overpressures. 

The theoretical “safe” altitude, is also given in table 1 2. A discussion of this parameter is 
given in appendix D. The safe altitude is defined to be the minimum altitude above the ground 
where cutoff can occur so that the disturbances which reach the ground wUl be acoustic in nature 
with significantly reduced overpressure intensity. The safe altitude is a strong function of the 
prev aili ng meteorological conditions. During the BREN tower tests it ranged from 400 m to about 
2200 m (1300 to 7200 ft) above the ground. 

Figure 22 shows the comparison between the airplane ground speed, Vq, and the shock 
propagation speed, Vp* . In theory, cutoff should occur where these two are equal. For cutoff 
above the ground, the airplane ground speed must be less than The data indicate this to be 

true for these rumble cases. As the airplane ground speed approaches Vp* , the moderate-to- 
heavy rumble and light boom cases predominated. For - Vq) > 6.0 m/sec (>19.7 ft/sec), 

low rumbles occurred. 

Figure 23 shows the comparison between the observed shock propagation speed determined 

from arrival times at the ground microphones, Vp^^, and the propagation speed as determined 

from the meteorological conditions, Vp* . The two agree very well, which indicates that these 

* rtiHx 

disturbances are propagating at the local speed of sound, Vp*, rather than at the slower airplane 
ground speed, Vq (see fig. 22). The range of (Vp^j^ - Vp*^^^^) is from -3.1 m/sec to +4.7 m/sec 
(-10.2 to +15.4 ft/sec), with all but five cases within ±2.0 m/sec (±6.6 ft/sec). 

The variation of the overpressure with cutoff altitude was also determined. This type of 
analysis would be required to determine an empirical ground speed “safety factor” for commercial 
threshold Mach number operations to avoid objectionable noise on the ground. Figure 24 shows the 
variation of the average overpressure of all 16 ground microphones, with the speed increment 
(Vp* - Vq). a groimd speed increment or “safety factor” of about 6 m/sec (19.7 ft/sec) is 
indicated if low rumbles are to be permitted, while a safety of over 10 m/sec (33 ft/sec) would be 
necessary for these test conditions if it were desired to produce no noise at the ground. 

It is interesting that some of the low rumble cases for (Vp*^^^^ - Vq) between 5 and 10 m/sec 
(16 and 33 ft /sec) have a higher average overpressure at the ground than the heavy rumble and low 
boom cases. This is most likely due to differences in rise times of the observed pressure signatures, 
with the more subjectively objectionable noise associated with the shorter rise times. Typical 
pressure signatures for these cases are sketched on figure 24. 


A theoretical safe altitude, Z§ was computed for each rumble-producing flight and was 
correlated with the subjective observations. The theoretical safe altitude computed is strictly valid 
only in the case of constant lapse rate of Vp*, so only order of magnitude comparisons could be 
expected. In general, the lapse rates observed during the tests varied considerably with altitude 
above the ground. Safe altitudes are given in table 12 for each flight. Figure 25 shows the variation 
of the overpressure on the ground with the difference between the theoretical cutoff altitude and 
safe altitude, (Z^ - Z^). For positive (Z^ - Z^) values, cutoff is higher than the safe altitude and low 
rumbles would be expected at the ground. When cutoff occurred below the theoretical safe altitude, 
heavy rumbles or low booms would be expected. Figure 25 shows that the transition to low rumbles 
for these data occurs at a value of (Z^ - Zg) of about +120 m (+400 ft). Thus, the safe altitude was 
underestimated for these data. The correlation indicated by these results is encouraging and 
indicates that the concept of a safe altitude may be of operational value. 


Ckound Reflection Coefficient 

Reflection of shock waves from the ground influences the observed overpressure magnitude. 
For regular reflection of weak oblique shock waves, is close to 2.0 (i.e., the observed 
overpressure is twice the incident overpressure). At the cutoff condition there is no reflection, so 
Kj^ should be near 1 .0. The manner in which varies close to cutoff was studied by Thomas 
(ref. 37) and Thery (ref. 38). These authors used the one-dimensional flow relationships to show 
that irregular reflection may occur close to cutoff where Kj^ may increase to 3.0 before decreasing 
to 1 .0 at cutoff. 

The method described in an earlier section was used to estimate Kj^ for the threshold Mach 
number flights. The resulting data were given in tables 9 and 1 1 of earlier sections. These data for 
categories Bl, B2, and B3 have been plotted in figure 26. In general, there is a large scatter in the 
results, which is to be expected since the method used to determine from the measurements is 
somewhat crude. The values for cases well before cutoff vary from the 0.95 to 2.2. There does 
not appear to be any increase in near cutoff, but rather a gradual decrease is suggested, 
beginning at about 0.175 rad (10^ from cutoff. This result suggests that the one-dimensional 
analysis by Thomas (ref. 37) and Thery (ref. 38) is not realistic. A two-dimensional analysis taking 
into account the axisymmetric nature of the actual flow may give better agreement with this 
experimental data. For the cases categorized as rumbles or low booms (table 12) varied from 
0.63 to 1.88, with an average value of 1.27. For some of these cases, however, the angle of 
incidence, 0, was appreciable due to an inversion in the Vp profile. This condition existed on 
October 23 and October 30, and Kj^ averaged 1.59 with a range of 1.41 to 1.72 for these six cases. 
The remaining rumble cases then average 1.20. 
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Analysis of Pressure Signatures 


The pressure signatures observed in the vicinity of caustics produced by threshold Mach 
number flight are significantly different from pressure signatures produced at higher Mach numbers. 
Nonlinear effects predominate at caustics, and the linear theory is not valid. Analysis of these 
caustic pressure signatures may be helpful in isolating these important nonlinear effects. 

Figure 27 shows the six “best” caustic pressure signatures observed, together with the linear 
theory steady, level pressure signatures at a slightly higher Mach number = 28.7 N/m^ 

(0.59 Ib/ft^)). Compared to the theoretical signatures these caustic signatures are of stronger 
intensity, longer duration, and exhibit the U-shape characteristic of caustics (rather than the 
well-known N-wave shape). The increase in overpressure for these six cases from the linear theory 
overpressure away from caustics ranges from 29% to 76%. It appears that small-scale turbulence 
may have produced spikes on the caustic pressure signature measured for pass 063 which resulted in 
an overpressure increase of 478%. Suspicion of this type of interaction stems from the sharpness of 
the spike (i.e., very small time duration) and the time delay between the initial pressure rise and the 
spike. Signatures of this nature are common for the temperature and wind conditions present near 
the ground during this series of flights. 

The duration of the theoretical pressure signature away from caustics is about 0.10 sec. The 
duration of the observed caustic pressure signatures, however, ranges from about 0.135 to 
0.140 sec, an increase of about 40% from linear theory. This increase in duration of the caustic 
pressure signatures can be caused by two effects— the aging caused by the long propagation path for 
the cutoff rays and nonlinear effects that become important very close to the caustic. 


Summary of Threshold Mach Number Data 

The analysis of the threshold Mach number data has provided valuable definitive information 
on the nature of shock waves produced by flight at speeds slightly greater than the threshold value, 
at the threshold value, and below the threshold value. At speeds greater than the threshold value, 
regular sonic booms are observed and linear sonic boom theory agrees very well with observed data. 
During several of the flights at speeds well above the threshold Mach number, however, caustics 
were observed. These caustics apparently were not produced by the cutoff condition but rather by 
relatively low-magnitude changes in the airplane speed. The maximum observed caustic intensity for 
these cases was 82.4 N/m^ (1.72 Ib/ft^), which gives an amplification factor of about 3 (compared 
to steady, level flight at slightly higher Mach numbers). 


Several flights produced caustics on the tower associated with the cutoff condition. These 
cases have helped to define the caustic phenomena more clearly. Significant increases in 
overpressure occur within about 100 m (330 ft) above and below the caustic. However, for these 
flights the caustic intensities were only about 25 to 50 N/m^ (0.5 to 1.1 Ib/ft^). This is an 
amplification of 1 to 1 .8 compared to steady, level flight. For one case, however, an overpressure of 
135.02 N/m^ (2.82 Ib/ft^) was measured where very sharp peaks characteristic of signature 
distortions due to microscale turbulence in the pressure signature occurred. Observed pressure 
signatures in the vicinity of caustics have a duration that is about 40% greater than pressure 
signatures produced at higher Mach numbers. 

Almost half of the threshold Mach number flights produced rumbles or low booms at the 
ground, since cutoff occurred above the tower. Analysis of the propagation speed of these 
disturbances at the ground showed that they propagate at the local sound speed, which was faster 
than the airplane ground speed. When low rumbles were produced on the ground, the airplane 
ground speed was at least 6 m/sec (19.7 ft/sec) lower than the maximum shock propagation speed. 
Comparison of a theoretical “safe altitude” for sonic boom cutoff (for which no objectionable noise 
would reach the ground) with the observed data was good considering the assumptions made in 
deriving it. 

The analysis of the ground reflection coefficient indicates a gradual decrease in Kj^ from 2.0 to 
about 1 .0 beginning at a shock wave angle of about 0.175 rad (10^ with respect to the vertical. 


ANALYSIS OF LONGITUDINAL 


ACCELERATION FLIGHT TEST DATA 


This section contains a discussion and presentation of the results of the analysis of the sonic 
boom measurements obtained from the 19 longitudinal acceleration flights. A brief description is 
given of the generation of the acceleration caustic, followed by presentation of the test 
measurements and their comparison with linear sonic boom theory. 


Introduction 

The caustic produced by a longitudinal acceleration will be produced by SSTs during the 
acceleration from subsonic to supersonic flight. In terms of sonic boom, a caustic is a cusp or fold 
that exists in a shock surface. During longitudinal accelerations through Mach 1 , a cusp is formed 
because the Mach number increases with time along the flight path. Figure 28 illustrates the 
development of shock waves and the caustic formation during acceleration from subsonic to 
supersonic Mach numbers. A single magnified sonic boom signature is observed at the caustic. As 
the acceleration continues, the caustic eventually reaches the ground, and is the first part of the 
shock wave to reach the ground. 

The caustic locus is the line which defines the progressive intersection of the caustic with the 
ground (see fig. 28) as it forms at progressively greater lateral distances from the flight path. In the 
region immediately surrounding the caustic locus, two or more pressure signatures are observed, one 
of which is somewhat stronger than the steady-flight signature but has nearly the same shape. At the 
caustic these two signatures are merged, forming a U-shaped signature. As the acceleration 
progresses, the caustic is reflected from the ground and the leading shock travels progressively 
farther ahead of the trailing shock until finally the trailing shock is completely refracted, so that 
only one pressure signature is observed. 

Two other U.S. flight test programs designed to provide measurements near caustics produced 
by longitudinal accelerations have been conducted prior to the 1970 BREN tower tests. The first 
tests were done in September and October 1961 at Edwards AFB (ref. 17). A total of seven 
longitudinal acceleration flights were performed; four of these cases are discussed in more detail in 
reference 16. A second test series consisting of five acceleration flights were conducted at Edwards 
AFB during October 1964 (ref. 18). These early tests, while of good quality, did not have the 
required microphone density to obtain definitive measurements near caustics. 


“Operation Jericho-Virage,” conducted by the French, provided the first definitive measure- 
ments of caustics. The resulting test data are documented in reference 26. These tests showed the 
necessity for using a dense network of microphones to observe the acceleration caustic. In addition, 
the theory of geometric acoustics was shown to be adequate for predicting caustic occurrences on 
the ground. 

Table 13 summarizes the 1 9 longitudinal acceleration flights conducted during the 1970 BREN 
tower tests. The nominal flight conditions for these flights were: 

Airplane altitude— 10.3 km (33 700 ft) 

Initial Mach number— 0.95 

Final Mach number— 1 .3 

Airplane heading— 035“ (true North) 

Acceleration magnitude— about 1.1 to 1.6 m/sec^ (3.7 to 5.1 ft/sec^) 

Table 13 contains the airplane weight, the maximum observed overpressure, boom time at 
microphone T-1, and an indication of the type of boom observed for each flight. About half the 
flights placed the acceleration caustic on or near the microphone array. 

In the next section the results of the comparison between theoretical and observed caustic 
locations are given. Then the sonic boom measurements are presented in three sections, according to 
nearness of the tower to the caustic-ground intersection. In the last section, sonic boom intensities 
predicted by the linear theory are compared with observed intensities near the caustic, and the 
nature of the caustic pressure signature is discussed. 


Caustic-Ground Intersection— Theory and Experiment 

A relatively simple method for calculating caustic locations numerically is described in 
appendix B of reference 33. Appendix C of reference 33 contains a comparison of the predicted 
caustic locations with observed caustic locations for five longitudinal acceleration flights from the 
French “Jericho-Virage” tests (ref. 26). The error of the predicted caustic locations ranged from 
-1 .37 to +0.69 km (-0.85 to +0.43 st mi). The error quoted in reference 26 for the French computer 
program “STROUM” for the same flights ranged from -2.0 to +2.4 km (-1.24 to +1.49 st mi). Both 
the method described in reference 33 and computer program “STROUM” are based on the linear 
theory of geometric acoustics, but the numerical calculation procedures are considerably different. 

In the calculation of theoretical caustic locations for the 19 BREN tower acceleration flights, 
the method of references 7 and 30 was used. This required the calculation of ray trajectories at 
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TABLE 13.-SUMMARY OF LONGITUDINAL ACCELERATION FLIGHT TEST CONDITIONS^ 
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small flight time increments. These theoretical caustic locations were then referenced with respect 
to the tower. Figure 29 shows calculated caustic locations for pass 045. Three calculated caustic 
locations are shown. One was calculated using rawinsonde data alone (S), another using all available 
meteorological data including the BREN tower and ATMS data (T), and the third using virtual 
temperature in place of ambient temperature for the T data (T^). The three different meteorological 
conditions change the caustic location slightly, but it should be remembered that T and T^ differ 
from S in the lowest 2 km (1.24 st mi) only. Thus, it would appear the upper level meteorological 
conditions are more important, particularly the wind conditions near the airplane altitude. In this 
case the wind was from the SSE, so that the ray directly beneath the airplane (0=0) was displaced 
slightly to the NW. 

Figure 30 shows the calculated and observed caustic locations for pass 093. For this case there 
was very good agreement between theory and experiment. The wind was generally from a northerly 
direction, so that the 0 = 0 ray was displaced to the SE. 

Table 1 4 gives the theoretical and observed caustic locations with respect to the tower for the 
remainder of these flights. The flights have been separated into categories according to nearness to 
the caustic. For seven flights the caustic passed over the tower and intersected the ground 
downtrack, so that only rumbles or low booms were observed at the tower. For three flights the 
caustic occurred well before the tower, so that regular sonic booms were observed. For eight of the 
flights, however, the acceleration caustic occurred on or near the microphone array. For these cases 
it is possible to make a more accurate comparison of theoretical and observed caustic locations. The 
errors in the theoretical caustic locations are summarized in table 14. The error ranges from 
-1.6 to +0.5 km (-1.0 to +0.3 st mi). A histogram of the errors is given in figure 31. The errors for 
passes 044, 045, 046, 047, and 049 are all negative, suggesting that the rawinsonde wind data near 
the airplane altitude are not quite accurate (the same rawinsonde observation was used for all of 
these cases). For these flight conditions a decrease of only 3.7 m/sec (12 ft/sec) in the tailwind 
component is required to cause the caustic to occur 1.0 km (3300 ft) farther downtrack. From the 
rawinsonde observation the tailwind component at the airplane altitude is about 18m/sec 
(60 ft/sec), so that a slight change in wind direction and/or magnitude could have produced such a 
change. 

The error in the theoretical caustic locations (for the seven cases when the caustic occurred on 
or near the microphone array) ranged from -1 .6 to +0.6 km (-1 .0 to +0.4 st mi). The range of error 
is very similar to the error for the five acceleration caustics of the Jericho-Virage tests referenced 
earlier (-1.4 to +0.7 km, -0.9 to +0.4 st mi). 
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TABLE 14.-C0MPARIS0N OF THEORETICAL GROUND CAUSTIC LOCATIONS WITH 

OBSERVED SHOCK WAVE CHARACTERISTICS-LONGITUDINAL ACCELERATION 
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For passes 041 and 043 no airplane radar data are available, and for pass 041 there are no sonic boom data. 


















































Comparison of Observed and Theoretical Boom Arrival Times 


For cases when sonic booms were observed on the tower, a comparison between theoretical 
and observed boom arrival times can be made. Theoretical arrival times can be calculated since the 
airplane position with time is known and since the sonic boom propagation time can be calculated 
for the ray trajectory which intercepts the tower base. The accuracy of the airplane position with 
time (i.e., radar accuracy), when translated to boom arrival time, is about ±0.2 sec. Errors of this 
magnitude can be attributed to the radar tolerance in the airplane position. Table 15 shows the 
comparison between observed and theoretical boom times. For two cases the arrival times at 
microphone G-14 were used, since the caustic (which is the first part of the sonic boom to occur on 
the ground) occurred past the tower and thus no theoretical arrival time at the tower base was 
predicted. The maximum error in the predicted arrival times for the 10 cases shown in table 15 is 
± 1 .0 sec. It may be significant that for all flights on August 28 the theoretical arrival times were 
too early, while for all flights on October 27 the theoretical arrival times were too late. This is 
consistent with the predicted caustic location errors shown in table 14, When an early arrival time 
was predicted, an “early” caustic location was also predicted (farther uptrack). The wind on 
August 28 was a tailwind, while on October 28 a headwind prevailed at the airplane altitude. 


TABLE 15.-COMPARISON OF OBSERVED AND THEORETICAL BOOM ARRIVAL TIMES 


Date 

Pass 

Bongo- 

pass 

Observed boom 
arrival time, 
PDT 

(mic. T-1) 

Theoretical boom 
arrival time 
at tower 
base, PDT 

Error, 

sec 

8*28 

040 

1-3 

0857:28.7 

0857:27.7 

-1.0 


045 

3-2 

1043:48.7 

1043:48.6 

-0.1 


046 

3-3 

1055:28.1 

1055:27.9 

-0,2 


047 

3-4 

1106:19.1 

1106:18.6 

-0.5 

10-27 

086 

1-1 

0843:47.2 

0843:47.7 

+0.5 


087 

1-2 

0853:31.9 

0853:32.7 

+0.8 


088 

1-3 

0904:13.7 

0904:14.3 

+0.6 


092 

3-1 

1146:53.6 

1146:54.4 

+0.8 


093 

3-2 

^1 158:36.4 
(G14) 

^1158:36.6 

(G-14) 

+0.2 

• 

094 

3-3 

^1210:51.0 
(G-14-Hl22m 
(400 ft) 

^1210:52.0 
(G-14-M22m 
(400 ft) 

+1.0 


^For these cases the caustic occurred downtrack of tower so arrival times 
at microphone G-14 were used. 
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Measurements Uptrack from Caustic-Ground Intersection 


Seven of the 19 acceleration flights produced caustics which passed over the tower, so that 
only rumbles and acoustic disturbances were observed (see table 13). For several of these cases it 
was possible to compute profiles of the acoustic fronts from measurements obtained on the BREN 
tower. Figure 32 shows the acoustic front locations (indicated by dashed lines) at the given 
reference times for passes 039, 044, and 049. The arrival time at microphone T-1 was used as the 
reference. In all cases the onset of the waves was used as a reference point to determine arrival 
times. 

The predicted and extrapolated caustic paths are indicated in figure 32. The caustic paths were 
estimated from the approximate caustic-ground intersections downtrack of the tower. Figure 32 
shows that the acoustic waves were oriented like a shock wave, and reflections from the ground are 
evident. The observed acoustic waves appear to be extensions of the caustic, which has a finite limit 
but produces acoustic disturbances that extend several thousand feet below it. Since the caustic is a 
moving source of disturbances, the acoustic wave would be oriented like a shock wave. For these 
cases it propagated over the ground (parallel to it) at a speed about 7% faster 20 m/sec (65.6 ft/sec) 
than the local sound speed near the ground. 

Figure 33 shows the pressure signatures observed on the tower for pass 049. While the 
signature looks very much like a typical N-wave, its time duration is about 20% greater than normal 
and it decays rapidly to an acoustic-type wave near midtower. 


Measurements Near Caustic-Ground Intersection 

Eight of the accelerations produced caustics on or near the microphone array. These 
measurements represent unique and valuable information on the nature of the caustic phenomena. 
Shock wave profiles are presented in this section for all of these cases, and observed pressure 
signatures on the tower and ground arrays are presented and discussed for most of them. 

Figure 34 shows the shock wave profiles for passes 043 and 092. For both of these cases the 
caustic intercepted the ground slightly uptrack of the ground microphone array. The observed 
pressure signatures on the ground during pass 043 are given in figure 35, while the tower data are 
given in figure 36. The caustic intersected the ground uptrack of microphone G-1, and the leading 
and trailing shock waves separated as the shock wave system propagated from microphone G-1 to 
G-1 4. On the tower (fig. 36) the incident leading shock is well defined, while the incident trailing 
shock is weak. The reflected shocks merge just above the tower top to produce the caustic. Some 
portions of pressure signatures have been omitted in figure 36 to fit them on a single page. 


Figure 37 shows the tower pressure signatures for pass 092. The incident trailing shock is 
somewhat better defined in this case and shows the typical U-shape characteristic of shock waves 
that have passed through the caustic. A moderate-intensity caustic can be seen at microphones T-9 
and T-10 where the reflected shocks have merged. 

Five of the caustic cases (passes 045, 046, 047, 093, and 094) are basically similar in character, 
since in each case the caustic was incident on the lower half of the tower. Shock wave profiles and 
observed pressure signatures are given in figures 38 through 46. In all of these cases the ground 
pressure signatures show acoustic-like disturbances occurring uptrack from the caustic-ground 
intersection. In most cases there is a rapid increase in overpressure within a few hundred feet of the 
caustic. Downtrack from the caustic the overpressure decreases with distance as the leading and 
trailing shock waves progressively separate. 

The tower pressure signatures and shock wave profiles for these five cases reveal a number of 
significant features. Above the caustic location on the tower the leading and trailing shocks are 
separated by a relatively small distance, which progressively decreases until they are merged at the 
caustic. Below the caustic, significant disturbances can occur for several hundred feet, but generally 
within 150 to 300 m (500 to 1000 ft) below the caustic (or after reflection from the ground), these 
disturbances are of relatively low magnitude and acoustic-like in nature. The pressure signatures 
near caustics exhibit the typical U-shape noted in previous flight tests. There is considerable 
variation in the maximum overpressures, however, and in signature shapes. 

Pass 046 is of particular interest since it produced the most intense caustic measured. The 
maximum overpressure at microphone G-8 was 240.4 N/m^ (5.02 Ib/ft^), and at microphone T-2 a 
maximum overpressure of 134.3 N/m^ (2.805 Ib/ft^) was recorded (for a ground reflection 
coefficient of 2.0 this becomes 268.6 N/m^ (5.61 Ib/ft^)). A distinguishing feature of pass 046 is 
that the distance between the leading and trailing shocks just above the caustic is significantly 
greater than for the other cases. The separation distance between the leading and trailing shocks 
increases with increasing airplane acceleration magnitude. This suggests that the acceleration may 
have been greater for pass 046, which resulted in the stronger caustic. 

To determine the effect of airplane acceleration on the caustic strength, airplane acceleration 
magnitudes were determined from the unsmoothed radar observations. A 7-sec flight time interval 
over which the shock waves near the caustic were generated by the airplane was used to calculate 
the average instantaneous acceleration. The distance between the leading and trailing shock waves is 
also related to the airplane acceleration. An arbitrary distance of 300 m (1000 ft) above the 
approximate caustic altitude was selected for measuring the shock wave separation, which was read 
directly from the shock wave profiles and is plotted in the upper half of figure 47. The maximum 
observed tower overpressures and shock separation distances are tabulated in table 1 6 and plotted in 
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figure 47. The data indicate a trend toward increasing caustic intensity with increasing acceleration 
magnitude. 


TABLE 16.-EFFECT OF ACCELERATION ON CAUSTIC STRENGTH 


Date 

Bongo- 

pass 

Pass 

Maximi 

obser 

jm overpressure 
ved on tower 

Shock 
separation 
at 300 m,(984 ft) 

Airplane 
acceleration 
from radar. A 

N/m^ 

Ib/ft^ 

Mic no. 

m 

ft 

m/sec^ 

ft/sec^ 

8-28 

3-2 

045 

109.6 

2.29 

T-2 

25 

82 

1.27 

4.2 


3-3 

046 

134.5 

2.81 

12 

92 

300 

1.48 

4.9 


3-4 

047 

93.4 

1.95 

T-4 

35 

115 

®1.25 

4.1 

10-27 

1-3 

088 

112.0 

2.34 

T-2 

50 

164 

1.34 

4.4 


3-2 

093 

107.3 

2.24 

T-7 

60 

197 

®1.36 

4.5 


3-3 

094 

80.0 

1.67 

T-5 

20 

66 

1.26 

4.1 


®For these flights the detailed radar data are suspect over the flight path interval of interest- 
therefore a longer averaging period was used. ' 


The last case (pass 088) is particularly interesting since it is anomalous when compared with 
other acceleration caustic cases. The shock profile for pass 088 is shown in figure 48 and indicates 
the nature of the anomaly. It appears that a caustic was produced uptrack of the tower (similar to 
pass 043, fig. 34). A second caustic, however, and leading-trailing shock wave system, also occur 
where in previous cases a trailing shock wave alone occurs. The measured pressure signatures for this 
flight are shown in figure 49. 

The presentation of the observed pressure signatures on the tower and ground, along with the 
calculated shock wave profiles, has shown that the caustic is formed at the cusp in the shock front 
where the leading and trailing shocks merge. The observed pressure signatures exhibit the 
characteristic U-shape, with sharp peaks at the bow and tail shocks. The maximum observed 
intensity was 240.4 N/m^ (5.02 Ib/ft^). There is some evidence to suggest that acceleration 
magnitude may affect the caustic intensity. 


Measurements Downtrack from Caustic-Ground Intasection 

For three of the acceleration flights (passes 040, 086, and 087) the caustic occurred well 
before the tower (about 2 to 5 km ( 1 .2 to 3. 1 st mi) from tower) so that a leading shock wave was 
observed, with a weak trailing shock wave following well behind the first shock wave. Figure 50 
shows the shock wave geometry system for pass 040 constructed from the tower arrival times. The 
leading and trailing shocks were separated by about 150 m (500 ft) at the ground. For this reference 
time the caustic occurred on the ground uptrack and passed above the tower, since it had been 
reflected earlier from the ground. The estimated caustic path is indicated. 
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Figure 51 is a reproduction of the pressure signatures observed on the ground for pass 040. 
The occurrence of “spikes” on many of the leading pressure signatures was believed to be due to the 
effects of small-scale atmospheric turbulence rather than airplane acceleration (refs. 20-22). The 
trailing signatures are much lower in magnitude and exhibit the Ushape characteristic of pressure 
signatures near caustics or those that have passed through it. The trailing signatures have thus been 
influenced by nonlinear effects. 

Figures 52 and 54 illustrate the shock wave profiles for passes 086 and 087, respectively. As in 
pass 040, the trailing shock wave is well behind the leading shock wave, with the caustic occurring 
well above the tower after reflection from the ground. Observed ground pressure signatures are 
given in figures 53 and 55. 

For these three passes it is possible to compute theoretical pressure signatures for the leading 
shock waves. The linear theory is not valid in the case of the trailing shock waves, since these have 
passed through the caustic and have been influenced by nonlinear effects. Figure 56 gives a 
comparison between calculated pressure signatures (Kj^ = 1 .0) and observed tower pressure 
signatures for passes 040, 086, and 087. In each case two tower signatures are given. The 
comparison between theory and experiment is reasonably good. In general, the calculated signature 
lengths and intensities agree fairly well. Considerable variation in the observed pressure signatures 
can be seen. These variations are believed to be caused by the effects of small-scale atmospheric 
turbulence not accounted for in the linear theory. The small-scale turbulence causes spikes and 
rounded signatures to occur alternately over small distances (refs. 20-22). 


Overpressure Variation Near Caustics— Theory and Experiment 

The linear theory of sonic boom propagation is not valid near caustics since, in the concept of 
ray-tube area from geometric acoustics theory, higher order effects have been ignored. The ray-tube 
area, A, is related to the sonic boom overpressure, AP, in the following manner: 

AP oc (8) 

Thus, near caustics where the ray-tube area becomes small, the overpressure becomes correspond- 
ingly large. At a caustic, where it is zero, an infinitely large overpressure is predicted, and the linear 
concept of ray-tube area is no longer vaUd. A fundamental consideration, then, is how close to 
caustics the linear theory is valid. A comparison with the experimental data is made in this section 
to empirically determine this limit. 
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Passes 045, 046, and 092 were selected for comparison with linear theory since for these cases 
caustics were produced before or near the BREN tower. Thus, overpressure data were available just 
downtrack of the caustic where the theory could be compared. Figures 57 through 59 give the 
overpressure variations near caustics for these cases. Two maximum overpressures are given, one for 
the leading shock, APj, and one for the trailing, caustic-like shock wave, AP 2 . At the caustic these 
are superimposed. The linear theory applies only to the leading shock wave. The data have been 
referenced to the caustic location on the ground for convenience. The accuracy of this location, 
however, is about 60 m (200 ft). Figures 57 through 59 also give a typical linear theory variation of 
overpressure close to the caustic. In addition, a steady, level flight reference intensity is given for 
comparison. For a ground reflection coefficient of 2.0 this reference intensity is about 56.5 N/m^ 
(1.18 Ib/ft^). Ground microphone numbers are indicated for the observed data. 

The overpressure variation for passes 045 and 046 (figs. 57 and 58, respectively) indicate that 
at about 400 m (1300 ft) from the caustic the nonlinear effects (AP 2 ) are more important than the 
linear effects (APj). That is, at 400 m (1300 ft) from the caustic APj is equal to or greater than 
AP 2 and thus the trailing, caustic-like signature begins to dominate. For pass 046 the linear theory 
predicts AP| fairly accurately within 100 m (330 ft) of the caustic. During pass 092 (fig. 59) a 
caustic was produced on the ground just uptrack of microphone G-1, and overpressure 
measurements farther downtrack of the caustic were obtained. Overpressure levels comparable to 
those produced during steady, level flight occurred at 400 to 500 m (1300 to 1600 ft) from the 
caustic. These data indicate that the theory gives reasonable results to within a horizontal distance 
of about 400 m (1300 ft) downtrack of the caustic-ground intersection. An examination of the 
tower pressure signatures (for example, fig. 42) shows that at a vertical distance of about 1 50 m 
(500 ft) above the caustic the linear and nonlinear effects are about equal. At distances closer to the 
caustic than these (within 400 m (1300 ft) horizontally downtrack of the caustic-ground 
intersection and within 150 m (500 ft) vertically) the nonlinear effects are important and the linear 
theory alone cannot be expected to give realistic results. These distances are somewhat dependent 
on the acceleration magnitude, generally increasing as the airplane acceleration magnitude decreases. 

In the analysis of the Jeiicho-Virage acceleration caustic data (ref. 26), a caustic “amplification 
factor” is given for the maximum observed amplification compared to a reference intensity 
measured away from the caustic. For the five longitudinal accelerations of the Jericho-Virage tests, 
this coefficient of amplification was found to be about 5. During the BREN tower acceleration tests 
the maximum observed overpressure was 240.4 N/m^ (5.02 Ib/ft^) at microphone G-8 on pass 046. 
The observed maximum overpressure at microphone T-2 for that same case, however, was 
134.3 N/m^ (2.805 Ib/ft^). For a reflection coefficient of 2.0 this becomes 268.6 N/m^ 
(5.61 Ib/ft^). The maximum caustic amplification observed, then, is this maximum overpressure 
divided by the reference steady, level overpressure value of 56.5 N/m^ (1.18 Ib/ft^). The resulting 
caustic amplification factor is 4.75, which agrees very well with the results of the Jericho-Virage 
acceleration tests. 
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Ray-Tube Area Limit Near Caustics— Linear Theory 


As shown in the previous section, at a distance of about 400 m (1300 ft) from the caustic the 
nonlinear effects are of the same order of magnitude as the linear effects. At this distance from the 
caustic the linear theory should still be valid, since the nonlinear effects do not predominate. Figure 
60 gives the variation of the ray-tube area with maximum overpressure and distance downtrack 
from the caustic-ground intersection. From figure 60 the ray-tube area at 400 m (1300 ft) is about 
2000 m (7000 ft) and the linear theory maximum overpressure is about 91.0 N/m^ (1.91b/ft^). 
Figure 61 compares the linear theory pressure signature 400 m (1300 ft) from the caustic with 
observed pressure signatures at that distance from the caustic. The linear theory predicts the leading 
shock wave intensity fairly well. Based on these data, a limiting value of the ray-tube area, A(S), of 
about 2000 m (7000 ft) is indicated. The area A(S) is the area normal to the ray, and is significantly 
different from the horizontal projection of the ray-tube area, Aj^, used in the basic theory. The two 
are related by the expression: 


n/ (9) 

sat^sin ©)Aj^ for small Uj^ 

Theoretical calculations at lower ray-tube area values than 2000 m (7000 ft) will give unrealistic 
results. 

By placing an empirically-determined limit on the ray-tube area it is possible to obtain realistic 
theoretical intensities at caustics. From figure 60 and using the caustic amplification coefficient of 
5, the ray-tube area limit would be about 300 m (1000 ft). A linear theory pressure signature was 
calculated using this ray-tube area limit. The maximum overpressure of this pressure signature is, as 
expected, 237.8 N/m^ (5.05 Ib/ft^). This linear theory “caustic” pressure signature is given in figure 
62, along with two observed caustic signatures for comparison. By using the ray-tube area limit the 
intensities compare favorably. The linear theory signature shape, however, does not agree since no 
nonlinear effects have been considered. This comparison shows that the nonlinear effects at caustics 
must be considered if realistic theoretical pressure signatures are desired. 


Discussion of Observed Phenomena Near Caustics 

A fundamental consideration is the nature of the nonlinear effects that transform the leading, 
incoming N-wave pressure signature into the U-shaped trailing pressure signature. The distance from 
the caustic where this process occurs is also of interest. In an effort to determine these the measured 
pressure signatures near caustics were analyzed in detail by subtracting out the theoretical N-wave 
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part of the measured signatures, leaving the nonlinear, caustic-like portion. When the leading and 
trailing shocks first begin to interact, they are superimposed in a linear manner. At the caustic, 
however, the superposition is nonlinear, since the N-wave is transformed to a U-wave. By 
decomposing the measured signatures into N-wave and U-wave components it may be possible to 
determine the nature and importance of the nonlinear effects with distance from the caustic. 

Figures 63 and 64 show typical results of this analysis for passes 045 and 046, respectively. In 
each case several signatures are analyzed at various distances from the caustic. The observed pressure 
signatures, the linear theory N-wave, and the nonlinear residual after subtraction of the N-wave 
from the observed signature are given for each signature. The N-wave is the calculated theoretical 
signature at about 400 m (1300 ft) from the caustic (A(S) = 2000 m (7000 ft)). 

The analysis in figure 63 for pass 045 shows the results for three signatures measured on the 
ground at about 50, 175, and 370 m (165, 575, and 1200 ft) downtrack from the caustic-ground 
intersection and for one tower signature. Qose to the caustic, the nonlinear residual consists of 
highly spiked shocks, with a nearly constant but slightly negative overpressure between the two 
shocks. The relatively large negative overpressure just prior to the tail shock is caused by the fact 
that the N-wave is of shorter duration than the observed signatures. The trailing U-wave generally 
has about a 10% to 20% longer duration than the incident N-wave. 

Figure 64 shows the analysis of the tower-observed signatures measured above the caustic 
during pass 046. These measurements in the vertical plane clearly show the merging of the leading 
N-wave and trailing U-wave. The observed pressure signature at microphone T-13 is well above the 
caustic where the N-wave and U-wave are separated and not interacting. At closer distances to the 
caustic the N-wave and U-wave are superimposed. As in the previous case, the nonlinear residuals 
exhibit the highly spiked shocks, and in all cases the nonlinear residual is well-behaved and 
reasonable. From this analysis it must be concluded that the N-wave and U-wave are superimposed 
in a linear manner up to very close to the caustic (within a vertical distance of about 30 m (100 ft)). 
Thus, the transformation of the incoming N-wave to the U-wave must also occur within this 
distance. This analysis method may prove to be useful for future investigators. 


Summary of Longitudinal Acceleration Data 

The analysis of the longitudinal acceleration data has shown that the method of linear 
geometric acoustics is capable of predicting caustic-ground intersections to an accuracy of about 
±1 .0 km (±3300 ft). This accuracy is about the same as that obtained for the Jericho-Virage data, 
which suggests that inaccuracies in the upper level meteorological conditions are responsible for 
most of the error. A comparison of observed and theoretical shock wave arrival times at the tower 


showed that they agreed to within ±1.0 sec. Most of this error is probably also due to lack of 
detailed knowledge of the meteorological conditions along the boom propagation path. These 
results indicate that linear geometric acoustics can be used with good accuracy when the airplane 
flight path and meteorological conditions are known. 

The analysis of shock wave profiles, airplane acceleration magnitude, and caustic intensity has 
shown that the caustic intensity may increase with increasing acceleration magnitude. Shock wave 
profiles are helpful in visualizing the phenomena, and it is clear that the caustic occurs where the 
leading and trailing shock waves merge, forming a cusp. Significant disturbances can occur below 
the actual caustic, since low-magnitude acoustic disturbances propagate several thousand feet below 
the caustic. 

The analysis of pressure signature intensities in the vicinity of caustics has shown that 
amplified intensities occur within about 400 m (1300 ft) downtrack of the caustic-ground 
intersection and within 150 m (500 ft) vertically above the caustic. Nonlinear effects become 
progressively more important as the caustic is approached. The maximum observed caustic 
overpressure for these accelerations was about 264 N/m^ (5.6 Ib/ft^). This is about five times the 
intensity that would be produced during comparable level flight at constant Mach number. Sonic 
boom pressure signatures calculated with linear sonic boom theory give reasonable agreement with 
the observations, but become unrealistic at distances closer than about 150 m (500 ft) vertically 
from the caustic. The separation of the observed pressure signatures near the acceleration caustic 
into linear and nonlinear components has shown that the transformation from the N-wave to the 
U-wave must have occurred within a vertical distance of about 30 m (100 ft) from the caustic for 
the case analyzed. 
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ANALYSIS OF LATERAL CUTOFF 


FLIGHT TEST DATA 


In this section the results of the analysis of the data for the nine lateral cutoff flights are 
discussed. Comparisons with theoretical results are shown for several cases. 


Introduction 

Lateral cutoff is similar to the threshold Mach number phenomenon; in each case the shock 
wave is completely refracted by atmospheric temperature and wind gradients. Lateral cutoff, 
however, occurs at all supersonic Mach numbers at the lateral extremity of the sonic boom carpet. 
The top part of figure 65 illustrates the lateral cutoff phenomenon. 

During the last decade there have been several experimental programs where sonic booms were 
observed at the side of the flight path (refs. 18-24). These pressure signature measurements have 
indicated that near the lateral extremities of the boom carpet, large variations in overpressure occur. 
These large variations are caused by atmospheric inhomogeneities that distort the shape of the 
signature as it travels long distances near the ground where atmospheric turbulence is most 
prevalent. Although large variations in overpressure occur near lateral cutoff, a general trend toward 
lower overpressure levels has been observed at the outer edges of the sonic boom carpet. 

The 1970 BREN tower sonic boom flight tests included nine flights designed to determine 
more exactly the characteristics of sonic boom near lateral cutoff. An illustration of the test 
arrangement is shown in figure 65. For the nine flights the nominal airplane flight conditions were: 
Mach number 1.3; altitude lO.OkmMSL (32 700 ft); and heading 125® true. The airplane lateral 
displacement from the BREN tower varied from 18.4 to 24.2 km (1 1.4 to 15.0 st mi). All nine of 
the passes were made on August 27, 1970 over a 3-hr period. Three flights were made, each 
consisting of three passes about 10 minutes apart. Table 17 summarizes the airplane flight 
conditions and observed sonic boom characteristics. The airplane weight was determined from 
pilot-observed fuel on-board at the “steady” and “breakoff’ points. The boom time at tower 
microphone T-1 is also given in table 17, along with the maximum observed overpressure recorded 
by any of the tower or ground microphones. 

The meteorological conditions during these lateral cutoff tests were important, since a fairly 
strong sidewind was present. At the airplane altitude the wind was from the SSW at a speed of 
about 15.2 m/sec (50 ft/sec), almost directly perpendicular to the flight path. This sidewind 
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produces an asymmetric shock pattern on the ground with respect to the flight path, as shown in 
figure 65 ; the shock-ground pattern is shifted to the NNW. The BREN tower was on the downwind 
side and thus in an area where lateral cutoff was extended due to the sidewind. 

Another important meteorological effect is the location of the maximum shock propagation 
speed, Vp. Normally this occurs at the ground. Occasionally, however, the maximum value occurs 
above the ground, so that an “inversion” exists in the Vp profile. When this is the case, lateral 
cutoff will occur at or above the maximum Vp altitude, which is above the ground. This condition 
can exist when there is a strong temperature inversion and/or when moderately strong tailwinds 
occur near the ground. During several of the BREN tower lateral cutoff flights an inversion did exist 
which somewhat complicated the data analysis. 


Lateral Cutoff Location— Theory and Experiment 

The actual location of lateral cutoff is somewhat difficult to determine from measured data. 
However, an excellent measure of how close the shock wave at the tower was to cutoff is the 
observed shock inclination angle relative to the vertical. Before cutoff a shock wave will be observed 
with a positive angle, 0, while after cutoff an acoustic front will be observed propagating faster at 
the lower levels (due to the normally higher temperatures and sound speed at the ground) so that 
the inclination angle, 6, is negative. Thus, 0 is a direct measure of closeness to lateral cutoff, and is 
useful for correlating measured data with theoretical calculations. 


The computation of the theoretical lateral cutoff location, was performed as follows: 

• The computer method of reference 30 was used to calculate lateral cutoff to within about 
0.0175 rad (1 .0“) of actual cutoff 


• The additional displacement to actual lateral cutoff (0= 0) was calculated by using the 
approximation (ref. 7): 


. 0(deg) 
57.3 


a(Z) 


• For cases when an inversion in the Vp profile existed, lateral cutoff was defined to be the 
location where 0=0 occurred anywhere on the tower. 


• To reduce inaccuracies caused by the boom-time-only meteorological observations at the 
tower, the values calculated for the three passes of each flight were averaged to 

obtain a more representative 


59 



Table 18 contains the theoretical lateral cutoff locations, the airplane lateral displacement 

from the tower, D, the observed shock inclination angle, 0, the ground reflection coefficient, Kj^, 
and the average ground microphones for each flight. 

Figure 66 gives an indication of the accuracy of the theoretical lateral cutoff location, 
compared to nearness to actual cutoff,0. The parameter (D-Yj^^^^) is plotted versus 0. is 

the theoretical lateral cutoff location with respect to the BREN tower. The agreement between 
theory and experiment is fairly good, and these calculations are accurate to about ±1.0 km 
(±3300 ft). It does not appear possible to achieve better accuracy than indicated in figure 66, since 
the meteorological conditions along the ray paths associated with boom at the tower are not known 
precisely. Thus, the accuracy indicated in figure 66 is probably an upper limit for these kinds of 
calculations. 

Another measure of the nearness of lateral cutoff is the shock propagation speed, Vpg (parallel 
to the flight path as determined from the propagation speed of the shock front normd to itselO 
compared to the airplane ground speed, Vq. In theory, lateral cutoff should occur when these two 
are equal. By analysis of the shock arrival times on the two-dimensional ground microphone array, 
it is possible to determine both of these parameters. Since microphones G-15 and G-16 were parallel 
to the airplane ground track, the shock arrival times, t^, can be used to obtain the airplane ground 
speed directly, as: 


^ (G-16)-tQ (G-15)] 

^mic AX 

where AX= 1081 m (3545 ft) = the distance between microphones G-15 and G-16. 

To obtain the propagation speed normal to the shock wave requires that the angle of the 
direction of shock propagation with respect to the line of microphones G-15 to G-16 be known. 
This can be calculated from: 

Vr . 

^mic 

[t^(G-14)-t^(G-l)]/AX 


where AX = 975.4 m (3200 ft) = the distance between microphones CM and G-14. 
The propagation speed, V_, is then: 


where Vp = a - Uj^; Uj^ = u cos (v-t)); 
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TABLE 18.-SUMMARY OF LATERAL CUTOFF FLIGHT TEST DATA ANALYSIS, AUGUST 27, 1970 
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‘Lateral cutoff occurred at midtower; therefore© is negative 
and estimated to be -0.034 rad (-2.0®). 



Table 1 9 presents these data and figure 67 is a plot of Vn • versus Vnc for six of the lateral cutoff 
. , '^mic 

flights (for three flights the disturbances were too weak to determine shock arrival times with 
accuracy). Also shown on figure 67 is Vpg*, the normal shock propagation speed with the effect of 
the presence of water vapor taken into account (see appendix A). The accuracy shown by these 
calculations is amazingly good, particularly for the Vpg* data. The values of 0, which indicate 
nearness to lateral cutoff obtained by an independent analysis, are also shown. The above indicate 
that for shock waves near lateral cutoff shock propagation speed is affected by water vapor and that 
linear geometric acoustics can be used to accurately predict the location of cutoff. 


Ground Reflection Coefficient 

A parameter of interest near lateral cutoff is the ground reflection coefficient, Kj^. A doubling 
of shock wave intensity occurs for oblique shock waves incident on smooth surfaces, so that 
Kr = 2.0. Measurements of K|^ for most natural ground surfaces have been between 1.8 and 1.9. As 
discussed in a previous section, however, Kj^ should approach 1 .0 near lateral cutoff where the 
shock wave is almost perpendicular to the ground. To provide further insight into this variation near 
lateral cutoff, the experimental data were analyzed to determine Kj^. A direct measurement of K|^ 
was not made, so the observed intensities on the tower and on the nearest ground microphones were 
used. The procedure was previously defined in the “Data Analysis Methods” section. 

Figure 68 shows the variation of K]^ withG. Before and close to lateral cutoff the data indicate 
that a Kj^ of about 1 .2 is representative, while past cutoff varied from 0.86 to 1 .08. Pass 030 is 
anomalous because cutoff and caustic-like signatures occured near midtower, while lower intensity 
signatures occurred on the ground, giving a K|^ value less than 1.0. As with the threshold Mach 
number measurements, these data do not indicate a significant increase in reflection coefficient near 
lateral cutoff. 


Overpressure Variation 

The variation of overpressure near lateral cutoff was studied using both theoretical predictions 
and the measurements. The theoretical variation is shown in figure 69. Directly beneat the airplane 
an overpressure of 57.0 N/m^ (1.21 Ib/ft^) is predicted, gradually decreasing with increasing lateral 
distance, and then increasing rapidly just prior to cutoff. The linear theory incorrectly predicts an 
infinite intensity at cutoff due to ray focusing. 

The data are plotted in figure 70. The average and extreme maximum overpressure, ^P^ax’ 
the pressure signatures observed on the ground are given for all nine flights. The shock wave 
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inclination angle, 0, was used as a measure of nearness to lateral cutoff. For each pass (except 034 
and 037) microphones G-1 through G-7 and G-8 through G-14 have been separated to give an 
indication of the overpressure variation along the microphone array. Microphone G-1 was closest to 
the airplane. 

The theoretical variation of overpressure with 0 near lateral cutoff is also shown in figure 70 
for two cases. The data for pass 030 are representative of the normal theoretical variation where an 
inversion in the Vp profile does not exist. The minimum overpressure calculated was about 
43.1 N/m^ (0.90 Ib/ft^) at the location where 0 was about 0.175 rad (10®). The theoretical 
variation of overpressure with © for pass 029, however, also shown in figure 70, shows a continued 
decrease in overpressure to much lower © angles. This is due to the effect of the inversion in the Vp 
profile that existed for pass 029. In plotting the data for 029 the minimum ©was used rather than© 
at the ground, since this value is more indicative of nearness to cutoff. It occurred at about 1 20 m 
(400 ft) above the ground. The better agreement of linear theory with the experimental data for 
pass 029 is significant, since it implies that the method is capable of calculating results to within 
0.035 rad (2.0®) of lateral cutoff for some combinations of meteorological conditions. 

Figure 70 indicates an interesting variation of overpressure near lateral cutoff. Relatively 
stronger intensities were measured very near lateral cutoff, and the pressure signatures for these 
passes (033 and 035) are caustic-like in nature. However, the highest overpressure measured was 
about 24.0 N/m^ (0.5 Ib/ft^). The average intensity for passes 033 and 035 was about 15.8 N/m^ 
(0.33 Ib/ft^). These relatively low-magnitude disturbances were perceived as low booms or rumbles. 

These results are much more definitive than an early study of overpressure variation near 
lateral cut off for these data, given in figure 9 of reference 24. The detailed analysis of nearness to 
lateral cutoff has added significantly to the interpretation of the results. 


Observed Pressure Signatures and Shock Wave Profiles 

The observed pressure signatures and shock front shapes are presented in this section. The 
results are separated into discussions of measurements obtained before, near, and beyond lateral 
cutoff. The shock wave profiles and tower pressure signatures have been presented on the same 
plots to the same distance scale. This required a transformation of the pressure signatures from the 
time coordinate to a distance coordinate by using the airplane ground speed (see “Data Analysis 
Methods” section). The pressure signatures observed on the ground are presented as observed 
(overpressure versus time). 
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Measurements before lateral cutoff. —Otly one of the nine lateral cutoff flights (pass 029) was 
performed with the airplane positioned so that lateral cutoff occurred past the tower. Figure 71 
shows the shock wave orientation in the vertical plane, along with several representative pressure 
signatures observed on the tower. Although the disturbance is of low magnitude (about 12.0 N/m^, 
0.25 Ib/ft^), a reflected shock wave is evident and the incident shock wave angle with the vertical, 
0, is about 0.070 rad (4®). The incident shock wave is very weak and iU-defined, while the reflected 
shock wave near the tower top has a well-defined signature shape. The reflected shock wave has 
significant curvature, since it was reflected from the ground earher and well before cutoff at a larger 
angle,0. 

For this case it was possible to calculate a theoretical pressure signature. The calculated 
signature is given in figiue 71 ; it has a maximum overpressure 17.72 N/m^ (0.37 Ib/ft^). It was also 
possible to calculate a theoretical shock front shape, which is also indicated. 

The pressure signatures observed at the 16 ground microphones are given in figure 72. There is 
considerable variation in the signature shape and intensity observed at each microphone. In 
addition, acoustic distrubances occur for a significant duration after the initial signatures. 

Measurements near lateral cutoff -Ihxte of the flights (passes 030, 033, and 035) produced 
shock waves at the BREN tower that were very close to lateral cutoff and exhibited caustic-hke 
pressure signatures. Figures 73 through 78 show shock front profUes and observed pressure 
signatures on the tower and ground for these cases. 

The tower data for pass 033 in figure 73 show a weak incident pressure wave, which may be 
the remnants of the incident shock wave, followed by the reflected shock wave. Since the incident 
wave is almost vertical, it is very near lateral cutoff. The pressure signatures on the ground in figure 
74 show one basic signature, followed by acoustic-like pressure disturbances. In general, the 
disturbances are acoustic in nature since the signatures are rounded, without shocks. 

Pass 035 shown in figures 75 and 76 is similar to pass 033 except that the incident wave is even 
weaker and occurs near the top of the tower only (indicated by the dashed line in fig. 75). The 
reflected pressure signature is caustic-like and is almost vertical at the ground. Both the pressure 
signatures on the bottom half of the tower and on the ground are caustic-like, but rounded and of 
relatively low intensity. Ground microphones G-8, G-9, and G-10 show an extra intermediate peak 
in the observed pressure signatures. 

Pass 030 is another case for which cutoff was near the tower. Figure 77 shows the shock front 
to be perpendicular to the ground from the middle of the tower to the tower top. Several pressure 
pulses can be seen ahead of the main part of the signature, and the signature itself is caustic-like 



with sharp peaks. Further evidence that lateral cutoff was just uptrack of the tower can be seen in 
the ground pressure signatures, shown in figure 78. These signatures have mostly random acoustic 
pressure variations, particularly microphones G-10 through G-14. Again an extra intermediate peak 
was observed at microphones G-8 and G-9, which suggests some effect caused by the small buildings 
at the tower base. 

Because cutoff occurred near the middle of the tower, pass 030 is somewhat anomalous when 
compared with other cases. For example, the ground reflection coefficient as calculated from the 
tower and ground overpressure intensities is less than one, since the tower intensities are greater 
than those on the ground. Also, the value of 0 at the ground was not obvious. It should have been 
negative since cutoff has already occurred at the ground. Thus, a 0 of -0.035 rad (-2.0®) was used for 
this case, based on the acoustic disturbances measured near the tower base. 

Measurements past lateral cutoff . of the nine flights produced disturbances at the tower 
that were well past cutoff, so that acoustic signals were observed. To observers on the ground, these 
flights (passes 031, 032, 034, 036, and 037) produced signals that were generally perceived as 
rumbles. The maximum overpressure of all these cases was 1 1.5 N/m^ (0.24 Ib/ft^); generally the 
overpressures were about 4.8 N/m^ (0.10 Ib/ft^). 

The data for pass 032 are given in figures 79 and 80. The acoustic front shown in figure 79 
propagates fastest at the ground, giving a negative 0. From the tower pressure measurements, it is 
difficult to distinguish a pressure signature as such. The ground pressure measurements show a rapid 
decrease of overpressure from microphone G-1 to microphone G-5. This case is anomalous in that 
respect, although it illustrates the large variations that can occur near lateral cutoff. The pressure 
measurements of microphones G-6 through G-14 are representative of the measurements of the 
remainder of these passes. 


Summary of Lateral Cutoff Data 

The nine lateral cutoff flights during the 1970 BREN tower tests have provided some 
interesting and valuable information about sonic boom characteristics near lateral cutoff. The most 
significant observation was the caustic-like signatures produced by ray focusing very close to the 
lateral cutoff location. The measured overpressure increases, however, were only about 50% higher 
than those observed just prior to cutoff and were much lower than those experienced directly 
beneath the airplane. The maximum overpressure observed for all the passes was about 24.0 N/m^ 
(0.50 Ib/ft^). Near lateral cutoff the incident shock waves were found to be very weak, while the 
reflected shock waves recorded on the tower in all cases were stronger. The ground reflection 
coefficient prior to lateral cutoff appears to be about 1 .2, while past lateral cutoff it is nearer to 1 .0 
but showed large variations. 
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Theoretical calculations of the lateral cutoff location using linear geometric acoustics are 
accurate to about ±1 .0 km (±3300 ft) for these cases, even though the exact meteorological 
conditions along the ray path were not precisely known. For the one case where it was possible to 
calculate a theoretical pressure signature, the agreement with the observed signature was reasonably 
good (the theoretical was 17.7 N/m^ (0.37 Ib/ft^) compared with the observed value of 

about 12.0N/m^ (0.25 Ib/ft^). Figure 81 is a summary of the shock front profiles and typical 
pressure signatures. 


CONCLUSIONS 


The 1970 sonic boom tests conducted at Jackass Flats, Nevada, have helped to define more 
exactly sonic boom phenomena associated with (1) flight near the threshold Mach number, (2) 
transonic acceleration, and (3) lateral cutoff. The use of the BREN tower was indispensable for 
measuring these test data, because it enabled vertical observations through the shock waves. Hence, 
it was not necessary to place caustics precisely at a point on the ground to obtain a dense sampling 
of measured pressure signatures. These vertical measurements also allowed the calculation of shock 
wave profiles and angles of incidence, which added substantially to the overall value of the data and 
to the interpretation of the test results. 

The observations near caustics are of particular value, since no methods are currently available 
for predicting realistic pressure signatures at caustics, because of the nonlinear effects that 
predominate there. Overpressure increases measured at caustics produced during threshold Mach 
number flight were relatively low, ranging from about 1.0 to 1.8 times those predicted for flight at a 
slightly higher Mach number. In one case where it appeared that small-scale atmospheric turbulence 
produced “spikes” on a caustic signature, the increase was higher. There is some evidence to suggest 
that the caustic intensity during threshold Mach number flight is affected by the shock propagation 
speed gradient. Caustics were also produced by small, inadvertent changes in the airplane speed 
during several of the “steady,” level threshold Mach number flights. These caustics were slightly 
stronger, with a maximum ampUfication of about 3. Measured overpressures at caustics produced by 
airplane accelerations ranged from 2 to 5 times those which would be observed during steady, level 
flight at about Mach 1.2. For these caustics, acceleration magnitude appeared to affect the 
intensity, with the strongest produced by the greatest acceleration. This conclusion is also 
supported by the data measured for some threshold Mach number flights, where small magnitude 
accelerations produced lower intensity caustics. 

Pressure signatures were also observed near lateral cutoff which resembled those measured at 
caustics. These disturbances were of very low intensity, however -less than one-half the intensity 
beneath the flight path. The meteorological conditions during the lateral cutoff flights were 
significant, since a relatively strong sidewind component prevailed and in several cases an inversion 
in the shock propagation speed profile existed at the tower. In theory, these conditions would be 
expected to have a strong effect on the measured overpressures. 

The maximum intensities and amplifications compared to the overpressures predicted directly 
beneath the airplane during steady, level flight are summarized in table 20. 
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TABLE 20-SUMMARY OF MAXIMUM INTENSITIES AND AMPLIFICATIONS 
PRODUCED BY VARIOUS FLIGHT CONDITIONS 


Type of caustic 

Maximum tower 
overpressure observed 

Approximate maximum 
amplification 


N/m^ 

Ib/ft2 

factor^ 

Steady Threshold 
Mach Number Flight 

,‘^50.8 

h.06 

1.8 

Low-Magnitude 
Acceleration Through 
Threshold Mach Number 

82.4 

1.72 

3 

Transonic Acceleration 

134.5 

2.81 

5 

Lateral cutoff 

^23.9 

'^0.5 

^1 


theoretical reference overpressure of 28.7 N/m^ (0.60 1b/ft^) was used as a 
reference for determining caustic amplification (steady flight, M =1.2, Kr = 1 -0). 

^For one case small-scale atmospheric turbulence appears to have produced 
“spil$£s" on a causjtic signature. This maximum overpressure was 135.0 
N/m^(2.82 Ib/ft^). 

^Since Kr 1.0, the amplification factor in terms of ground-obseryed 
boom is 0.5. 


Analysis of pressure signatures near caustics and comparison with linear theory showed that 
nonlinear effects predominate at caustics. Nonlinear effects become important within about 150 m 
(500 ft) vertically above the acceleration caustic (or about 400 m (1300 ft) horizontally from the 
caustic-ground intersection) and within 100 m (330 ft) above and below the threshold Mach 
number caustic. Caustic pressure signatures are quite different in shape from normal sonic booms, 
since they typically exhibit a U-shape rather than an N-wave shape and the intermediate shock wave 
(for the F-104) is generally not present. In addition, caustic pressure signatures have a significantly 
longer duration than normal N-wave pressure signatures. These measurements of sonic boom 
characteristics at caustics produced by several different flight conditions should prove to be of great 
value in extending theoretical methods to include calculations near caustics. 

During the threshold Mach number flights and the lateral cutoff flights, considerable 
information was obtained on the acoustic disturbances that occur past the cutoff condition. 
Generally, these disturbances propagate at the local speed of sound rather than the airplane ground 
speed. By detailed analysis of these data it was possible to isolate the effect of the presence of water 
vapor on the propagation speed. An increase was noted, in agreement with theory. For low 
supersonic flight and when the increase in propagation speed is 1.0 m/sec (3.3 ft/sec) or more, this 
effect should be taken into account. 




Additonal analysis of the rumble data produced during flight near the threshold Mach number 
showed that ‘‘low rumbles” occurred when the airplane ground speed was at least 6 m/sec 
(19.7 ft/sec) lower than the maximum shock propagation speed. This, roughly, is the “safety 
factor” (or reduction in allowable airplane ground speed to avoid objectionable noise at the ground) 
that would be required for these flight and meteorological conditions. Comparison of a theoretical 
“safe altitude” for sonic boom cutoff with the observed data was good considering the assumptions 
made in deriving it. 

A characteristic of sonic boom disturbances near the cutoff condition is the presence of 
“precursors” or pressure pulses that propagate ahead of the basic pressure signature. Precursors were 
produced during the threshold Mach number flights and the lateral cutoff flights when shock waves 
were near the cutoff condition, and thus propagating at speeds close to the local propagation speed 
parallel to the ground. 

The ground reflection coefficient was calculated for a number of cases, and the data indicate a 
gradual decrease from about 2.0 to 1 .0 near the cutoff condition. This is in direct contrast to several 
theoretical studies which predict an increase to 3.0 near cutoff. 

Comparison of theoretical calculations with the observed data showed good agreement in all 
cases where it was possible to perform such calculations. Shock wave intensities, inclination angles, 
and arrival times were calculated and compared with the observations. The accuracy of the 
theoretical calculations was limited by the accuracy of the input airplane flight path and 
meteorological data, particularly for the low supersonic flights where the shock wave travels long 
distances almost parallel to the ground. Inaccuracies in the upper level wind conditions are 
considered to be an important source of error. Caustic locations during transonic acceleration and 
lateral cutoff locations can be predicted to within ±1.0 km (3300 ft). Shock wave intensities agree 
reasonably well, as do signature shapes when the effects of small-scale turbulence are neglected in 
the observed data. Shock wave arrival times can be predicted to within ±1 .0 sec. Thus, the 
comparison of theory and experiment has tended to verify the theory and has also indicated its 
range of validity. 
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APPENDIX A 


EFFECT OF HUMIDITY ON SHOCK PROPAGATION SPEED 


The presence of water vapor has a second-order effect on the speed of sound and the shock 
wave propagation speed. Normally in sonic boom studies this effect can be neglected. For shock 
wave propagation near cutoff, however, where shock waves are deteriorating to acoustic waves, this 
effect may be important. The small change in velocity caused by humidity acting over relatively 
large propagation distances for cutoff rays can result in significant alteration in the direction of 
sound waves or low-magnitude sonic boom waves. 


The effect of humidity on sound velocity is attributed to the fact that the density of water 
vapor is approximately 60% of the density of dry air. Thus, humid air is slightly less dense and has a 
slightly higher speed of sound than dry air at the same temperature and pressure. The normal 
method for accounting for the effect of humidity is to use a “corrected” temperature called the 
virtual temperature. 


The virtual temperature, Ty of moist air is given by the expression: 

X =x 

'v ITTTr 

where: 


(Al) 


T = temperature of moist air 

r = mixing ratio of moist air (ratio of mass of water vapor to the mass of dry air with 
which the water vapor is associated) 

e = ratio of molecular weight of water vapor to that of dry air (0.62198) 

Ty = temperature which dry air must have at the given pressure in order to have the same 

density as moist air when the densities of both dry and moist air are deduced from 
the equation of state for ideal gas 

Tables for the calculation of Ty are given in reference 39 for various pressure, temperature, and 
humidity conditions. 



A useful form of equation (A1 ) is: 


D = (T,-T)=[4-i] [t^]t 


D" (0.60777) [fiTT)] T 

where T is expressed in *1C and D is expressed in “K or °C. 

Some of the results given in earlier sections have indicated that it may be possible to use the 
experimental data to determine the increase in sound speed because of the presence of water vapor. 
The lateral cutoff and threshold Mach number flight data for cases where acoustic propagation was 
occurring would appear to be useful for this purpose. Table 21 is a comparison of the shock 
propagation speed as determined from the arrival times at the ground microphones, Vp^.^ and the 
maximum shock propagation speed calculated from the meteorological conditions on the tower, 
V . Both of these calculations have been discussed earlier. The only cases considered were those 
for’^^ich it was possible to calculate a reasonable accurate propagation speed from identifiable 
pressure signature features. In addition, only cases for which acoustic signals were observed were 
used in this analysis. On two days (October 23 and 30), significant inversions existed in the Vp 
profile well above the tower, so that it was not possible to determine accurately a representative 
shock propagation speed. In most cases the maximum occurred within 30 m (100 ft) of the ground. 
Reference 34 contains Vp data for each pass. 

The difference between Vp . and Vp on the tower is a measure of the increase in shock 
propagation speed caused by water vapor. 9^he data in table 21 show considerable variation for 
individual passes. This is to be expected, since the accuracy of both Vp^.^ and Vp are the same 
order of magnitude as their difference. (The increase in shock propagation to be isolated is less than 
1.0 m/sec (3.3 ft/sec) or about 0.3% of the propagation velocity itself.) Daily averages of the data 
are also given in table 21 ; these indicate less variability since for only one day is (Vp^.^ - Vp^^^) 
negative. Combining all August measurements and all October measurements by a second averaging 
gives an increase of 1.1 m/sec (3.6 ft/sec) for the August data and 0.25 m/sec (0.8 ft/sec) for the 
October data. 

The theoretical increase in shock propagation speed for the observed atmospheric conditions is 
given in table 22. Average humidity and temperature conditions were used. An average increase 
during the August period was 1.0 m/sec (3.3 ft/sec) and during October 0.20 m/sec (0.7 ft/sec). The 
lower temperatures and humidity during October give a much tower increase in shock propagation 
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TABLE 21. -COMPARISON OF OBSERVED SHOCK PROPAGATION SPEED OVER GROUND WITH 
PROPAGATION SPEED FROM METEOROLOGICAL OBSERVATIONS 
(Acoustic Non-Inversion Cases Only®) 
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speed than during August. The experimental data averaged over August and October agree very 
closely with the predicted shock propagation speed increase. Figure 82 shows the comparison 
between theory and experiment. While considerable scatter is evident in the data, the averaged 
experimental data agree surprisingly well with theory. These data indicate that the observed change 
in shock propagation speed due to water vapor is approximately equal to the 

theoretical increase in shock propagation speed. 

In view of these results, it would be desirable to account for the increase in shock propagation 
speed due to water vapor during low supersonic Mach number flight when the water vapor content 
is high enough to cause an increase of more than 1.0 m/sec (3.3 ft/sec) in the shock propagation 
speed. For supersonic flight at speeds appreciably higher than the threshold Mach number and for 
low-temperature, low-humidity conditions the effect need not be considered. Figure 83 shows the 
theoretical increase with relative humidity for an atmospheric pressure condition of 1000 mb 
(2088.5 Ib/ft^) and several temperatures. For the propagation speed to change by more than 
1.0 m/sec (3.3 ft/sec), temperatures in excess of about 12.0“C (53.6’’F) are required with high 
relative humidity. At a temperature of 30.0" C (86.0" F) a relative humidity of only 35% is required 
since the air is capable of holding much more water vapor before saturation. 

Harris (ref. 40) has presented the results of some interesting laboratory work on the effects of 
humidity on the velocity of sound in air, obtained at two frequencies in the lower audio range. His 
results indicate that as moisture is added to dry air the velocity of sound at first decreases until the 
minimum velocity is reached at 14% relative humidity. Above about 30% relative humidity the 
velocity increases linearly with increasing moisture content. The initial decrease in sound speed with 
increasing relative humidity is attributed to the effect of molecular relaxation. As humidity is first 
added to dry air, the decrease in velocity because of molecular relaxation is larger than the increase 
in velocity because of density changes. His measurements were performed in air at a constant 
temperature of 20.00° C (68.0" F) and at normal atmospheric pressure. The estimated accuracy of 
the sound speed increase is ±0.05 m/sec (±0.16 ft/sec). 

Harris’ work raises an important question with regard to sonic boom waves propagating at 
near-sonic speeds. It is generally acknowledged that the humidity at which the molecular absorption 
is maximum (a decrease in sound speed) increases with increasing frequency. This suggests two 
things when sonic boom signals (consisting of a wide range of frequencies) are considered. First, the 
higher frequencies would tend to have slower propagation speeds because of higher molecular 
absorption, while the lower frequencies would tend to have higher speeds because of lower 
molecular absorption. Thus, the higher frequencies would tend to lag behind the lower frequencies, 
after propagation over large distances. Secondly, for low-frequency signals the molecular absorption 
is at a minimum, so that the full increase in propagation speed predicted by theory would be 
expected. 


Harris has also done some experimental work on the attenuation of sound (ref. 41). In this 
study he found that a 2000-Hz sound wave was attenuated to 0.3 its original ampUtude after 
traveling 400 m (1300 ft). A 300-Hz signal, however, is attenuated the same amount after traveling 
4 km (2.5 St mi). Thus, the higher frequency waves are attenuated much more rapidly with distance. 
This attenuation (and slowing) of higher frequency sound waves suggests that past sonic boom 
cutoff where acoustic disturbances occur, the higher frequency components would be lost, leading 
to rounded, smooth signatures. 

Hodgson and Johannesen (ref. 42) give an approximate expression for the thickness (rise time) 
of weak, fully dispersed shock waves. They found that humidity has a large effect on the shock 
thickness. For an overpressure of 24.0 N/m^ (0.5 Ib/ft^) the thickness is calculated to be between 
6 m (19.7 ft) (dry air) and 6 mm (0.24 in.) (hot humid conditions) at 300® K (80.3® F) and between 
3 m (9.8 ft) and 30 mm (1.18 in.) at 270® K (26.3® F). Thus, even before cutoff where shock waves 
exist, shock thickening may occm producing rounded pressure signatures (see fig. 71). 
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APPENDIX B 


AIRPLANE F-FUNCTION DATA 


To obtain accurate theoretical sonic boom pressure signatures, it is necessary to compute the 
airplane F-function. The F-function specifies the initial conditions for the sonic boom propagation 
and is calculated from the detailed geometry and lift distribution of the airplane. The F-function 
was derived by Whitham (ref. 43), and references 8, 44, and 45 describe methods for computing it. 
To facilitate the computation of the effects of airplane Mach number, load factor, weight, and the 
ray azimuth angle, the format developed by NASA Langley Research Center was used, where the 
geometry and lift components are nondimensionalized. In general, the F-function varies most 
significantly with airplane station, L, Mach number, M, ray azimuth angle, 0j., and lift coefficient, 
Cj^. It is represented as: 


F = Ff(FA + FBCOS0f) 


(Bl) 


= F,Fi 


where Fj- is that part that varies significantly with Mach number, F^ is the volume contribution, 
and Fg is the lift contribution. Mathematically they are: 


Ff=M2/(2/3)'/2 


(B2) 


/(J/2)(nLW/q)\ 

(LpiW— ) 


The parameters F^j and Fgj are the nondimensionalized volume and lift contributions, 
respectively. The other parameters are: 

Lp = a basic length used to nondimensionalize the F-function 

M = Mach number 

n^ = lift load factor 

Pq = atmospheric pressure at the airplane altitude 


q = 1/2 Pq ~ dynamic pressure 

W = airplane weight 

/3 = (M^- 1)1/2 

y = 1.4 = ratio of specific heats for air 

0J. = ray azimuth angle 

Tables 23 and 24 give the nondimensionalized geometry and lift components of the 
F-function, F^j and Fgj, for Mach numbers of 1.1 and 1.3 for the F-104 airplane. Since the 
parameters F^| and Fgj vary slowly with Mach number (most of the Mach number dependence is 
contained in the Ff factor), they may be used over a fairly wide range of Mach numbers. In the 
computer method (refs. 7 and 30) used to calculate the sonic boom pressure signatures the 
parameter F^ is calculated for each flight time and Mach number. The Mach 1.1 F-function data 
were used for the threshold Mach number pressure signature calculations and the Mach 1.3 
F-function data were used for the acceleration and lateral cutoff pressure signature calculations. 

The data in tables 23 and 24 are given at various airplane stations, L/Lp, beginning at the nose 
and extending past the tail. The F-function reference length is 0.3048 m (1.0 ft). The units of F^^ 
and Fg are (ft)^/^. The airplane reference wing area, Sg^pp, was 18.708 m^ (196.0 ft^). The 
airplane weight, W, determines the lift contribution to the F-function. 

Typical F-function data are given in figures 84 and 85 for the particular flight conditions noted 
on the figures. The contributions of various airplane components to the F-function are indicated in 
figure 85. 
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Table 23.-f-function data for f-io4 at mach i.io 


^Lp = 0.3048 m (1.0 ft); S^pp = 18.208 (196.0 ft^) 


L.ft 

>=A1 


L,ft 

Fai 

Fbi 

0.0 

0.0 

0.0 

36.0 

-0.314922 

0.0083620 

0.5 

0.0325166 

0.0 

37.5 

- 0.243227 

0.0060987 

1.0 

0.0568285 

0.0 

38.5 

- 0.091745 

0.0039997 

1.5 

0.0753496 

0.0 

39.0 

0.001296 

0.0020510 

2.0 

0.0895217 

0.0 

40.0 

0.2136 

- 0.0029230 

3.0 

0.1082570 

0.0 

40.5 

0.0872 

- 0.0058339 

3.5 

0.1139427 

0.0 

41.0 

0.0585 

- 0.0097013 

4.5 

0.1209745 

0.0 

41.5 

0.0522 

- 0.0151971 

5.5 

0.1303826 

0.0 

42.0 

0.0271 

- 0.0110773 

7.5 

0.1571810 

0.0 

42.5 

0.0086 

- 0.0078617 

8.5 

0.1731034 

0.0 

43.0 

- 0.0031 

- 0.0062546 

9.0 

0.1782168 

0.0 

43.5 

-0.0114 

- 0.0052313 

9.5 

0.1733862 

0.0 

44.0 

-0.01755 

- 0.004504 

10.0 

0.1586302 

0.0 

44.5 

- 0.02157 

- 0.003952 

10.5 

0.1350020 

0.0 

45.0 

- 0.00785 

- 0.003517 

11.5 

0.0672169 

0.0 

45.5 

- 0.00269 

- 0.003165 

12.0 

0.0256514 

0.0 

46.0 

- 0.00125 

“ 0.002873 

12.5 

- 0.0142418 

0.0 

46.5 

- 0.00187 

- 0.002627 

13.0 

- 0.0459985 

0.0 

47.0 

- 0.01911 

- 0.002417 

14.0 

- 0.090001 1 

0.0 

47.5 

- 0.00714 

- 0.002236 

15.0 

- 0.1143158 

0.0 

48.0 

- 0.01119 

- 0.002077 

17.0 

- 0.1402204 

0.0 

48.5 

- 0.01630 

- 0.001937 

18.0 

- 0.1520094 

0.0 

49.5 

- 0.02972 

- 0.001703 

18.5 

- 0.1575704 

0.0 

50.5 

- 0.04805 

- 0.001513 

19.0 

- 0.0976976 

0.0 

51.1 

0.14020 

- 0.001428 

19.5 

- 0.0103902 

0.0 

52.5 

0.05361 

- 0.001228 

20.5 

-0.0158538 

0.0 

53.15 

- 0.00702 

- 0.001156 

21.5 

- 0.0179638 

0.0 

54.0 

0.12815 

- 0.001070 

22.0 

0.0248128 

0.0 

54.5 

0.15800 

- 0.001025 

22.5 

0.0592704 

0.0 

54.9 

0.17490 

- 0.000996 

23.0 

0.0341000 

0.0 

55.5 

0.02278 

- 0.000944 

23.5 

0.0131156 

0.0 

56.0 

- 0.14983 

- 0.000907 

24.0 

- 0.0066519 

0.0 

56.5 

- 0.25893 

- 0.000867 

25.0 

- 0.0416093 

0.0 

56.9 

- 0.19683 

- 0.000850 

26.5 

- 0.0721385 

0.0 

57.15 

- 0.21905 

- 0.000833 

27.0 

- 0.0751240 

0.0 

58.0 

- 0.02270 

- 0.000782 

27.5 

- 0.0748659 

0.0 

58.5 

0.03222 

- 0.000755 

28.0 

- 0.0724169 

0.0 

59.0 

0.05335 

- 0.000730 

28.5 

- 0.0677570 

0.0030731 

59.45 

0.14208 

- 0.000706 

29.0 

0.1054490 

0.0052179 

60.0 

0.06701 

- 0.000684 

29.5 

0.2103205 

0.0066011 

60.5 

0.04495 

- 0.000663 

30.0 

0.2679040 

0.0076445 

61.5 

0.02953 

- 0.000623 

30.5 

0.2422265 

0.0084255 

62.5 

0.02153 

- 0.000588 

31.0 

0.1714775 

0.0090480 

63.5 

0.01753 

- 0.000555 

32.0 

0.0091100 

0.0124630 

64.5 

0.01485 

- 0.000526 

32,5 

- 0.0749998 

0.0131268 

65.0 

0.01391 

- 0.000512 

33.0 

-0.135454 

0.0132260 

, 67.5 

0.01003 

- 0.000452 

34.0 

-0.211430 

0.0123497 

70.0 

0.00719 

- 0.000402 

35.0 

- 0.269263 

0.0104246 

72.0 

0.00600 

- 0.000369 

35.5 

- 0.296427 

0.0092520 

75.0 

0.00426 

- 0.000327 


«For Lp = 1.0, Fa = Fai and Fg = Fg, [| ( 
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TABLE 24.-F-FUNCTION DATA FOR F-104 AT MACH 1.30 

^Lp = 0.3048 m (1.0ft);SRpp = 18.288 (196.0 ft^) 


L,ft 


«= A 1 


' B 1 


uft 


•= A 1 


^B1 


0.0 

0.0 

0.5 

0.025489 

1.0 

0.046537 

1.5 

0.063821 

2.0 

0.077941 

3.0 

0.098322 

4.0 

0.110278 

5.0 

0.120083 

6.0 

0.132273 

7.0 

0.146507 

8.0 

0.162150 

8.5 

0.189494 

9.0 

0.171036 

9.5 

0.163363 

10.0 

0.147299 

11.0 

0.093184 

12.0 

0.021022 

13.0 

- 0.035122 

14.0 

- 0.071911 

15.0 

“ 0.094154 

16.0 

- 0.110615 

17.0 

- 0.125486 

18.0 

- 0.138975 

18.5 

- 0.076983 

19.0 

- 0.028085 

20.0 

- 0.026429 

21.0 

- 0.025551 

21.5 

0.026107 

22.0 

0.040255 

22.5 

0.024784 

23.0 

0.009648 

23.5 

- 0.005522 

24.0 

- 0.020072 

24.5 

- 0.034020 

25.0 1 

- 0.046156 

26.0 

- 0.062328 

27.0 

- 0.067120 

28.32 

- 0.062800 

28.5 

- 0.060170 

29.0 

- 0.055389 

29.5 

0.022966 

30.0 

0.195671 

30.5 

0.304881 

31.0 

0.306149 

31.5 

0.232487 

32.0 

0.126461 

32.5 

0.026559 

33.0 

- 0.070748 

33.5 

- 0.154147 

34.0 

•- 0.2061 60 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0033474 

0.0056780 

0.0071764 

0.0082951 

0.0091224 

0.0099622 

0.0119982 

0.0136219 

0.0140847 

0.0139967 

0.0135073 

0.0126914 


35.0 

36.0 

36.5 

37.0 

37.5 

38.0 

38.5 

39.0 

39.5 

40.0 

40.5 

41.0 

41.5 

42.0 

42.5 

43.5 

44.5 

45.0 

46.0 

47.0 

48.0 

50.0 

50.5 

51.0 

51.5 

52.0 

52.5 

53.0 

53.5 

54.0 

55.0 

56.0 

57.0 

58.0 

59.0 

60.0 
61.0 
62.0 

63.0 

64.0 

65.0 

66.0 
68.0 

69.0 

69.5 

70.0 

71.0 

72.5 

74.0 

75.0 


- 0.282686 
- 0.350166 
- 0.352058 
- 0.329608 
- 0.283138 
- 0.202292 
~ 0.099395 
0.017605 
0.147210 
0.290607 
0.039094 
0.018795 
0.003757 
- 0.006566 
- 0.013909 
- 0.024085 
- 0.031073 

- 0.033937 

- 0.038558 

- 0.042389 

- 0.045358 
- 0.049865 

- 0.050753 
- 0.012266 

0.086256 

0.073575 

0.063988 

0.057613 

0.050119 

0.044532 

0.036273 

0.029772 

0.025192 

0.021930 

0.018686 

0.015523 

0.013558 

0.011948 

0.018012 

0.009658 

0.008575 

0.007879 

0.006461 

0.005985 

0.005630 

0.005535 

0.005053 

0.004599 

0.004105 

0.003798 


0.0102549 

0.0082707 

0.0074282 

0.0066129 

0.0058070 

0.0043933 

0.0022456 

- 0.000349 

- 0.003303 

- 0.006570 
- 0.011211 
- 0.017391 
- 0.010724 

- 0.007794 

- 0.006257 

- 0.004527 

- 0.003540 
- 0.003186 

- 0.002644 

- 0.002290 
- 0.001982 
- 0.001564 
- 0.001478 
- 0.001400 
- 0.001329 
- 0.001263 
- 0.001203 
- 0.001148 
- 0.001097 
- 0.001050 

- 0,000965 

- 0.000891 

- 0.000827 
- 0.000770 
- 0.000719 

- 0.000674 

- 0.000633 

- 0.000596 

- 0.000563 

- 0.000533 

- 0.000505 

- 0.000480 

- 0.000435 
- 0.000416 

- 0.000406 

- 0.000397 
- 0.000380 
- 0.000357 
- 0.000336 
- 0.000323 


a 


0 /"l w 



APPENDIX C 


THEORETICAL RELATIONSHIPS NEAR 
CUTOFF DURING THRESHOLD MACH NUMBER FLIGHT 


A number of useful relationships can be derived relating various parameters near the cutoff 
condition. For reference, several of these are presented in this section. 


Threshold Mach Number 


The threshold Mach number has been defined to be the maximum airplane Mach number for 
which complete shock wave refraction will occur above the ground. Since the refraction is caused 
by gradients of temperature and wind, the threshold Mach number is dependent on the sound and 
wind speeds at the airplane altitude and at the cutoff altitude. It can be derived from the equation 
describing the change of the wave normal direction with altitude. This relationship directly beneath 
the airplane for steady, level flight is given by equation ( 1 7) of reference 7 as: 


where: 


COS0 = 


- a(Z) 

Co + “n ® 


(Cl) 


0 = inclination angle of wave normal below the horizontal 

a(Z) = sound speed at altitude Z 

u,.|(Z)= wind component at altitude Z parallel to flight path (tailwind is negative) 


Co = (-^ -Uti ) = (Ma»-u„ ) = Vr 

o sinp^ o G 


(C2) 


The subscript o indicates initial values at the airplane altitude. 


The condition for cutoff or complete refraction is that 0 = 0 or cos 0= 1 .0. Then we have 

a(Z) 


cos0= 1.0 = 


[Mxao-Un^] + u^fZ) 


(C3) 
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Solving for M j gives 


Mt = 3^ 

The parameter |a(Z) - u 
between the airplane and the ground. 

Nicholls (ref. 46) has stated that the above equation for Mj (and the equivalent equation in 
ref. 32) is incorrect since no sidewind effect has been taken into account. In reality, however, 
Nicholls’ equation agrees with previous results when the differences between coordinate systems for 
the wind terms are considered. Nicholls failed to realize that the wind component Uj^ is the horizontal 
wind component in the vertical plane of the wave normal. Nicholls’ equation (12) of reference 46 
which was used for comparison between his and previous results is in error since Uj^ is not equal to 
his wind term u. The wind component Uj^ is defined as u cos For 0 = 0 this becomes 

u cos ( 0 - 1 ?) where: 

u = horizontal wind speed 

0 = airplane heading (true North) 

77 = direction from which wind blows (true North) 

u = heading angle of wave normal (true North) 


j^(Z)|max maximum shock propagation speed, that occurs 


Theoretical Relationships Near Cutoff During Threshold Mach Number Flight 

In deriving relationships for propagation near cutoff it is useful to start with the basic 
relationship describing the inclination of the wave normal below the horizontal, 0 , with sound and 
wind speed differences between the airplane and a particular altitude. This is given by equation 
(Cl). Two other parameters of interest are the airplane ground speed, Vq, and the shock wave 
propagation speed, Vp. These are given by 

= (C5) 

Vp = a(Z) - Un(Z) (C 6 ) 
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Also: 


n- 

® (M a„ - u„^) t u„(Z) - Vg + u„(Z) 

“T'S^Kmax^X} 

By algebraic manipulation of the above equations it can be shown that the following relationships 
exist as long as 0 is greater than 0®. This means that cutoff has not occurred. 


|Vp + u_(Z)| Vp 
” |Vq + Ujj(Z^ *Vq 

(C7) 


(C8) 

<VG-Vp)=a(Z)[<'^“*®'] 

(C9) 

<M-MT)-(VG-Vp_^^)/a„ 

(CIO) 

Vq = {a(Z>/cose}-Up(Z> 

(Cll) 


Other forms of the above relationships could be derived, but in this study these were found to be 
useful. For example, since cos 0 was available from the tower shock arrival times and a(Z) and 
Uj|(Z) were known from the tower meteorological observations, the airplane ground speed, Vq, 
could be calculated from equation (Cll). The value of (M - M-p) is a useful parameter, as is 
(Vq - Vp). They apply only as long as cutoff has not occurred, however. In general, for 0 greater 
than 0“the following will be true: 

cos0<l.O; (M - Mj)>0.0; Vp/VQ< 1.0. 
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APPENDIX D 


THEORETICAL SAFE ALTITUDE FOR SONIC BOOM CUTOFF 


Below the caustic produced by threshold Mach number flight the signal decreases in amplitude. 
Several years ago Dr. Wallace D. Hayes determined a theoretical distance below the cutoff altitude 
where the signal will be reduced to an acoustic disturbance of relatively small intensity. Dr. Hayes 
based his analysis on the theory of reference 10. Two parameters are required which include the 
signal length, L, and a relative radius of curvature, R, of the caustic (relative to the ray). The ratio 
L/R must be small. 

According to the theory the safe altitude, Zg, is 

Zg = (R)*/^ (L)2/3 (D1) 


where R depends on the meteorological conditions. 


R = 


a 

3(a-u„)/aZ 


(D2) 


It should be noted that in this analysis R should be approximately constant with altitude and 
always positive. R can be negative under certain atmospheric conditions such as in very stable lapse 
rate conditions or with a tailwind increasing with altitude (negative Uj^). In physical terms, a 
negative R means that cutoff cannot occur. 


It is of interest to calculate a typical safe altitude for the F-104 airplane used in the BREN 
tower tests. A typical pressure signature length is about 39.6 m (130.0 ft) at the ground for the 
flight test conditions. For standard no-wind conditions, R is calculated to be about 87 600 m 
(288 000 ft). Thus, the safe altitude is 517 m (1700 ft). For the atmospheric conditions during the 
BREN tower tests the safe altitude calculated from the local meteorological conditions using 
equation (Dl) ranged from about 400 m (1300 ft) to over 2000 m (6600 ft) above the ground. 
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FIGURE 1.-TEST SET UP DURING BREN TOWER TESTS 
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FIGURE 6-REFRACTION OF SONIC BOOM SHOCK WAVE FOR FLIGHT NEAR 
THE THRESHOLD MACH NUMBER 
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FIGURE 12.-SHOCK WAVE PROFILE AND TOWER PRESSURE SIGNATURES, 
' M> Mj-PASS 117, OCTOBER 30, 1-2 
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FIGURE 16.-SHOCK WAVE PROFILE AND TOWER PRESSURE SIGNATURES 
M - M-r-PASS 058, AUGUST 31, 2-4 
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FIGURE 17.-SHOCK WAVE PROFILE AND TOWER PRESSURE 
SIGNATURES, M Mj-PASS 063, AUGUST 31, 4-1 
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FIGURE 18. -SHOCK WAVE PROFILE AND TOWER PRESSURE SIGNATURES, 

PASS 102, OCTOBER 28, 3-2 
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FIGURE 21--VARIATI0N OF OVER PRESSURE NEAR CAUSTIC 
-THRESHOLD MACH NUMBER FLIGHT 
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SHOCK WAVE INCLINATION ANGLE, 0, DEG 

FIGURE 26. -VARIATION OF GROUND REFLECTION COEFFICIENT NEAR CUTOFF-THRESHOLD MACH NUMBER 
FLIGHT, M > M-r 


o 



115 











FIGURE 28-DEVELOPMENT OF SHOCK WAVES DURING ACCELERATION FROM SUBSONIC MACH NUMBER 
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FIGURE 29 -COMPARISON OF THEORETICAL AND OBSERVED CAUSTIC 
LOCATIONS-PASS 045 
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FIGURE 30.-COMPARISON OF THEORETICAL AND OBSERVED 
CAUSTIC LOCATIONS-PASS 093 
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TOWER TOP 

FIGURE 33.-TOWER PRESSURE SIGNATURES UPTRACK FROM CAUSTIC-GROUND INTERSECTION 
PASS 049, AUGUST 28, 4-2 
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IFIGURE 35.-GROUND PRESSURE SIGNATURES JUST DOWNTRACK OF CAUSTIC-GROUND 
INTERSECTION-PASS 043, AUGUST 28, 2-3 
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FIGURE 36.-TOWER PRESSURE SIGNATURES JUST DOWNTRACK OF CAUSTIC-GROUND 
INTERSECTION-PASS 043, AUGUST 28, 2-3 
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MAXIMUM OBSERVED OVERPRESSURE NEAR DISTANCE BETWEEN LEADING AND TRAILING SHOCKS 

CAUSTIC ON TOWER, APp, LB/ FT^ AT 990 FT ABOVE CAUSTIC, S, FT 
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FIGURE 51.- GROUND PRESSURE SIGNATURES DOWNTRACK FROM CAUSTIC-GROUND 
INTERSECTION- PASS 040, AUGUST 28, 1-3 
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FIGURE 56 -CALCULATED AND OBSERVED PRESSURE 
SIGNATURES-PASSES 040, 086, 087 
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FIGURE 62.-LINEAR THEORY PRESSURE SIGNATURE WITH RAY-TUBE AREA 
LIMIT COMPARED WITH CAUSTIC SIGNATURES 
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FIGURE 64.-SEPARATI0N OF OBSERVED PRESSURE SIGNATURES 
NEAR THE ACCELERATION CAUSTIC INTO LINEAR 
AND NONLINEAR COMPONENTS, PASS 046, AUGUST 28, 3-3 
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FIGURE 65.-EXPERIMENTAL SETUP FOR THE LATERAL CUTOFF FLIGHTS 
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FIGURE 66-COMPARISON OF OBSERVED LATERAL CUTOFF LOCATION WITH THEORY 
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, FIGURE 67.-NORMAL SHOCK PROPAGATION SPEED NEAR LATERAL CUTOFF 
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FIGURE 68.-VARIATION OF GROUND REFLECTION COEFFICIENT, Ko, NEAR LATERAL CUTOFF 
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FIGURE 69.-LINEAR THEORY VARIATION OF OVERPRESSURE WITH LATERAL DISTANCE 
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FIGURE 70 .-VARIATION OF MAXIMUM OVERPRESSURE ON GROUND NEAR LATERAL CUTOFF 
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FIGURE 71. -SHOCK WAVE PROFILE AND TOWER PRESSURE SIGNATURES 
BEFORE LATERAL CUTOFF-PASS 029, AUGUST 27, 1-1 
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FIGURE 73.-SHOCK WAVE PROFILE AND TOWER PRESSURE SIGNATURES NEAR LATERAL 
CUTOFF-PASS 033, AUGUST 27, 2-2 
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FIGURE 75.-SHOCK WAVE PROFILE AND TOWER PRESSURE SIGNATURES NEAR 
LATERAL CUTOFF-PASS 035, AUGUST 27, 3-1 
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FIGURE 77.-SHOCK WAVE PROFILE AND TOWER PRESSURE SIGNATURES 
NEAR LATERAL CUTOFF-PASS 030, AUGUST 27, 1-2 




1400r 


1200 b 


O 1000 

fsl 

UJ 

CO 

< 

CO 


800 


5 
o 

H 
ui 
> 

O 

OQ 

< 

LU 

O 

D 

t 600 

t: 

< 


400 


200 



-800 -600 -400 -200 0 200 400 

DISTANCE FROM TOWER, D, FT 

FIGURE 79.-ACOUSTIC WAVE PROFILE AND TOWER PRESSURE SIGNATURES PAST 
LATERAL CUTOFF-PASS 032, AUGUST 27. 2-1 
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FIGURE 81 .-SUMMARY OF SIGNATURES AND SHOCKWAVE PROFILES NEAR LATERAL CUTOFF 





THEORETICAL CHANGE IN PROPAGATION SPEED DUE TO WATER VAPOR, (Vp* • Vp), FT/SEC 

FIGURE 82.-CHANGE IN PROPAGATION SPEED DUE TO WATER VAPOR-THEORY 
AND EXPERIMENT 
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FIGURE 83.-INCREASE IN SHOCK PROPAGATION SPEED DUE TO WATER VAPOR-THEORY 
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FIGURE 84.-F-FUNCTION FOR F-104 AIRPLANE, M = 1.3 
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FIGURE 85.-F-FUNCTION FOR F-104 AIRPLANE, M = 1.10 
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